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SECTION  5 


ERROR  RATE  ESTIMATION  UTILIZING  CHANNEL  MEASUREMENTS 
BASED  UPON  RECEIVED  SIGNALS  ALONE 


This  section  of  the  report  deals  with  the  specification  and 
analysis  of  techniques  for  predicting  error  rates  for  the  fading 
channels  of  interest.  A variety  of  techniques  are  considered 
varying  in  complexity  and  performance.  It  has  been  pointed  out 
in  Section  2 that  the  IDS  channel  should  be  analyzed  on  a quasi- 
stationary basis  because  bad  performance  occurs  only  during  deep 
fades  which,  by  link  design,  are  (hopefully)  made  to  be  rare.  It 
is  more  meaningful,  for  channel  quality  monitoring  of  LOS  channels, 
to  estimate  the  error  rate  at  each  instant  of  time  due  to  the  ob- 
served channel  strength  and  frequency  selectivity,  than  to 
average  over  an  assumed  stationary  fading  channel  because  the 
fading  is  too  slow  to  obtain  meaningful  averages  over  the  time 
duration  of  the  fading  state.  Thus,  the  results  here  should  be 
regarded  as  applicable  to  the  HF,  Troposcatter , and  Satellite 
Scintillation  channels  although  strictly  for  comparison  purposes 
we  have  presented  average  error  rate  estimates  for  the  LOS  chan- 
nel also.  When  the  fading  is  flat  over  the  typical  DCS  LOS  band- 
width of  14  MHz,  a very  reliable  estimate  of  quasi-stationary 
error  rate  is  obtainable  by  using  the  measured  |T(f,t)|2  in  an 
appropriate  error  rate  formula.  The  existence  of  frequency  selec- 
tivity may  be  determined  from  measurements  of  |T(kF,t)|^  but 
quasi-stationary  error  rates  cannot  be  estimated  knowing  |T(f,t.)|^ 
alone.  Special  probing  signals  are  required  to  estimate  T(f,t). 

The  discussion  below  is  divided  into  six  sections.  Sections 
5.1 -5.3  analyze  three  classes  of  techniques  for  average  error  rate 
estimation  due  to  flat  fading  in  diversity-combining  receivers. 

An  appropriate  "instantaneous"  error  rate  is  determined  and  then 
averaged  over  the  fading  for  the  techniques  of  Section  5.1.  In 
Section  5.2,  the  probabilities  that  instantaneous  SNR's  are  below 
thresholds  are  measured  and  related  to  the  desired  error  probabi- 
lities. For  Section  5.3  a formula  is  assumed  relating  average 
error  rate  to  SNR  and  to  diversity  branch  correlation.  Measure- 
ments of  the  latter  are  used  in  the  formulas  to  estimate  error  rate. 
Section  5.4  presents  examples  of  the  estimation  of  irreducible 
error  rate  due  to  fast  fading  and  frequency-selective  fading. 
Section  5.5  r < nsiders  the  estimation  of  error  rate  for  an  FDM-FM 
data  transm: ..sion  system.  Section  5.6  discusses  the  error  rate 
estimation  problem  for  advanced  anti-multinath  modems  which  employ 
in-band  diversity.  Finally,  Section  5.7  deals  with  the  analysis 
of  error  rate  estimation  techniques  when  interference  is  present. 

In  all  analyses  complex  Gaussian  fading  statistics  are  assumed, 
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5*1  Flat  Fading  Error  Rate  Estimation  by  Averaging 
''instantaneous"  Error  Rates 


5.1.1  Introduction 

In  this  section  we  address  the  problem  of  estimating  the 
error  rate  of  a flat  fading  channel  by  averaging  instantaneous 
error  rate  estimates.  Two  techniques  are  examined.  The  first 
technique  uses  estimates  of  the  instantaneous  SNR  on  each  diver- 
sity branch  to  form  error  rate  estimates  for  each  of  the  diversity 
branches.  The  branch  error  rate  estimates  are  then  used  to 
estimate  the  channel  error  rate,  assuming  independently  fading 
diversity  channels.  The  second  technique  that  is  examined  is  to 
use  instantaneous  SNR  estimates  to  estimate  the  SNR  at  the  output 
of  the  diversity  combiner.  The  combiner  output  SNR  is  used  to 
estimate  an  error  rate  at  that  instant.  These  instantaneous 
error  rates  are  averaged  to  form  an  estimate  of  the  channel  error 
rate.  This  technique  does  not  require  the  assumption  of  indepen- 
dently fading  diversity  channels.  However,  it  is  more  difficult 
to  estimate  error  rates  when  the  modem  error  rate  is  low.  A 
technique  is  proposed  for  alleviating  this  problem  by  amplifying 
the  apparent  error  rate  and  extrapolating  to  the  correct  value. 

Finally,  a technique  is  studied  for  predicting  short-term 
error  bursts . 

5.1.2  Use  of  Diversity  Branch  Error  Rate  Estimates  for 
Error  Rate  Prediction 

In  this  section  we  consider  the  problem  of  estimating  the 
error  rate  of  a flat  fading  channel  from  estimates  of  the  instan- 
taneous error  rate  on  each  diversity  branch,  assuming  independently 
fading  diversity  branches.  Figure  5.1  is  a block  diagram  of  the 
system  to  be  analyzed.  The  estimation  procedure  is  as  follows. 

The  data  signal  is  transmitted  over  the  diversity  channel,  picked 
off  at  RF  (or  IF)  in  the  receiver,  and  sent  to  a magnitude-squared 
channel  transfer  function  estimator  of  the  type  analyzed  in 
Section  4.1.  The  output  of  the  estimator  is  sampled  and  combined 
with  an  estimate  of  the  noise  power  to  form  an  estimate  of  the 
instantaneous  SNR.  From  this  instantaneous  SNR  estimate,  we  ob- 
tain an  estimate  of  the  instantaneous  error  rate  for  that  diversity 
branch.  These  are  averaged  to  form  an  estimate  of  the  average 
error  rate,  which  are  combined  with  the  average  error  rate  esti- 
mates of  the  other  diversity  branches  to  give  a channel  error  rate 
estimate.  A possible  correction,  for  the  bias  in  the  error 
probability  estimator  is  shown  in  dashed  lines. 
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Figure  5.1  Block  Diagram  of  Error  Rate  Estimator 


For  analysis  purposes,  a 2-phase  DCPSK  system  is  used.  Thi3 
scheme  is  sufficiently  similar  to  other  signaling  schemes  so  that 
the  results  of  this  section  can  be  used  for  comparison  purposes. 
For  a 2-phase  DCPSK  system  and  Gaussian  noise,  the  error  rate 
conditioned  upon  the  SNR  is  given  by 


(5.1) 


where  y is  the  SNR.  From  (5.1)  and  estimates  of  the  signal-to- 
noise  ratio  on  the  diversity  branch  and  in  the  pc^  sampling 
instant,  we  can  form  an  estimate  of  what  the  ensemble  error  rate 
would  be  on  this  branch  and  at  this  sampling  instant  if  the  other 
branches  were  not  there: 


(5.2) 


A ^ til 

where  is  an  estimate  of  the  SNR  on  the  k branch  and  at  the 

pcft  instant.  Assuming  stationarity , the  error  rate  on  the  k*-*1 

branch  can  be  formed  by  averaging  the  above  estimates.  This  gives 

an  estimate  of  the  error  rate  for  the  kth  branch  acting  alone: 


A 


(5.3) 


For  many  diversity  combining  techniques,  it  is  possible  to 
express  the  error  rate  of  the  channel  approximately  in  terms  of 
the  error  rates  of  the  particular  branches  acting  alone  when 
independent  fading  occurs  on  the  diversity  branches.  For  example, 
with  M the  order  of  diversity,  we  can  express  the  error  rate  by 
[5.1] 


where 


M 

Pe“CfT  Pk 

k-1 


(5.4) 


C 


for  equal  gain  combining 

for  selection  combining 
for  maximal  ratio  combining 
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Equation  (5.4)  is  exactly  true  for  maximal  ratio  combining,  asymp- 
totically true  for  low  error  rates,  but  still  useful  for  error 
rates  up  to  10” 2 to  10“^  for  the  other  combining  techniques. 

From  (5.4)  we  can  estimate  the  error  rate  by 


M A 

K - c IT  \ 

k=l 


(5.5) 


The  ability  of  (5.5)  to  estimate  the  error  rate  will  be 
determined.  From  (4.250)  of  Section  4.3,  the  estimate  of  the 
SNR  on  the  k^h  diversity  branch  and  at  the  p^1  sampling  instant 
can  be  approximately  represented  by 


<1+V 


Y- 


kp 


+ 


(5.6) 


where  ykp  is  the  actual  instantaneous  SNR,  and  ekp  and  6kp  are 
independent  estimate  errors  with  first-  and  second-order  error 
moments  approximately  given  by 


€ 


kp 


0 


kp 


2 _ 2 _2 
kp  CS  "k 


(5.7) 


where,  for  the  estimator  of  Figure  4.1,  e and  are,  respectively, 
given  by  (4.57)  and  (4.100).  In  (5.7),  Tk  is  the  mean  SNR  of  the 
kth  diversity  branch.  Using  (5.6)  in  (5.3),  the  estimate  of  the 
error  rate  on  the  diversity  branch  acting  alone  is  given  by 


P,  - 


2K 


K 

L 

p=i 


■[(1+Vr 


kp+6kp] 


(5.8) 


Averaging  over  the  instantaneous  SNR  estimation  errors  and 
channel  fluctuations  gives 


(5.9) 
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Since  the  estimation  errors  due  to  the  channel  selectivity,  5kp, 
are  independent  of  the  errors  due  to  noise  and  data,  eko,  r 


(5.10) 


For  zero  mean  complex  Gaussian  channels,  has  a one-sided 
exponential  distribution  with  mean  Therefore,  performing  the 
average  over  the  channel  fluctuations  gives 


ir.  + l + e,  rT  , 

k kp  kj 


(5.11) 


The  channels  studied  have  allowed  sufficiently  accurate  esti- 
mation of  the  magnitude  squared  of  the  channel  transfer  function 
so  that  E'jPk}  can  accurately  expressed  in  terms  of  the  first- 
and  second-order  moments  of  the  error.  Therefore,  we  can  approxi- 
mate (5.11)  by 


E 


2 2_2  - 
f c 2 f '^k 

1+T  rk  + H 

(tv  + id 


(5.12) 


Noting  that  for  the  system  considered  we  have  [5.1]  for  the 
average  error  rate  of  the  kfch  diversity  branch  acting  alone 


, = ___i 

k 2(rk+l) 


(5.13) 


and  upper  bounding  slightly  the  term  in  by  assuming  F^/ (F^+l) = 1, 
it  immediately  follows  that 


2 

rk + 


e{pJ  - p 


(5.14) 


Hence,  the  diversity  branch  error  rate  estimator  is  biased.  The 
quantity  e is  known  and  o will  normally  be  measured.  If  e| 
could  be  estimated,  the  bias  could  be  corrected.  While  insuffi- 
cient  time  was  available  to  study  this  question,  it  appears  that 
tg  can  be  measured. 

The  error  variance  in  estimating  will  be  determined  next. 
Squaring  (5.8)  and  averaging  gives  for  the  second  moment  of  the 
estimate  of  the  error  rate  for  the  kt^1  diversity  branch  acting 
alone 


E 


- 


— E I EL't(1+fkP)ykP+6kP+(1+fkq)ykq+6kq]i 

4K2  p=l  q=l 


(5.15) 


To  evaluate  the  above  expectation,  the  joint  density  of  y, 
and  y,  must  be  known.  For  a complex  Gaussian  channel,  the  p 
joint  density  is  a function  of  the  correlation  coefficient  of 
the  channel  complex  envelope.  In  particular,  defining 


P 


pq 


A 

;Tg(F,tp)Tk(F>ta)} 

■H 

[!Tk(Ot)l2j 

[ 

(5.16) 


where  T^Cf,!)  is  the  time-varying  channel  transfer  function  of  the 
k*-b  divers  ity  branch,  then  the  joint  density  of 
expressed  as  [5.2] 


^kp  and  ^kq  can  be 


P^'kp’rkq) 


% 


m 


°\  ^ - 


ykpykq 


^kp^kq 

r.  (1-p2  X 

i_k  Mpqj 


(5.17) 


pq/ 


Using  (5.17)  to  perform  the  average  over  y,  and  y.  gives  for 
the  second  moment  of  the  error  rate  estimate  on  the^ke“  diversity 
branch  acting  alone 
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K K 


E{$  ' i £ S E 


"^6kp+6kq^ 


4K‘  p-1  q-1  |rEa+ckpXi+ckqXi-Ppq)  +rk(2+ckp+fkq)  + 1 


(5.18) 


Expressing  in  terms  of  the  first-  :ad  second-order 

moments  of  the  errors  gives 


e-Tp^V  - 


2 2 2 

1 + 2<srk  + f 


lkr  4K(2r.+l) 

2 2 2 

l + c^  + 2e  K K 

+ M= £ £ 


4 K 


P--1  <V  1)2  - (Ppqrk)2 


(5.19) 


where  the  errors  at  different  sampling  times  were  assumed  to  be 
independent. 

To  simplify  notation,  we  define  ^*(Ppq)  as 


* (ppqrk) 

*<PpJ  ■ ^ 2 

pq  <rk+1>  - w 


(5.20) 


where  we  note  that  <I,(Ppq)~>0  as  Ppq“*0.  Figure  5.2  presents  typi- 
cal time  dependencies  for  the  functions  discussed  in  this  section. 

Substituting  the  above  expression  for  <S(pPq)  into  Eq.  (5.19) 
and  using  the  expression  for  P,  given  by  (5.13;  gives  approxi- 
mately for  E|P|~j 


1 . O 2-2  , 2 

1 + 2gsrk+< 
4K(2rk  + l) 


(i^2  + 2e2) 

K 

K r 1 

\ S k ! 

4K2 

E 

p=l 

Eji-tPpqJ  + l] 

p^q 


(5.21) 


■tr^^ ^>wpww  _ txwwqr* i^ow 


Using  (5.12)  and  (5.21)  to 
error  rate  estimate  of  the  kth 


evaluate  the  variance  of  the 
diversity  branch  acting  alone  gives 


a . 


2 2 2 

1 + 2*sriS 

4K(2rk  + l) 


1 J.  2_2  , . 2 

l + gSrk  + 2g  2 
K Fk 


1 + esrk  + 2f2  o K K 
S~T P ^ E *(P 

p=l  q-1 

p^q 


pq' 


(5.22) 

whore  higher  than  second-order  error  moments  were  considered 
negligible. 

For  stationary  scatter  channels,  ppq  can  be  expressed  in 
terms  of  |tp“tqj  and,  for  equally  spaced  samples , pp(~  can  be  ex- 
pressed in  terms  of  |p-q| . Thus,  for  this  case,  thedouble  sum 
of  (5.22)  can  be  written  as  a single  sum  giving 


0 0 0 2 9 9 

1 + 2*Srk  + < 1 + *5rk  + 2*  .2 

4K(2rk  + l)  K k 


+ 


i , 2 2 2 

1 + fsrk  + 2* 

K2 


K-l 

E 


r=l 


(K  - r)  <S>[p(r)] 


(5.23) 


where  P(r)=nPpq  with  r=  |p-q|. 

When  K is  large,  $[p(r)]  contributes  to  the  above  sum  only 
when  K-r»K.  Furthermore,  the  second  term  is  negligible  for  all 
cases  of  practical  interest.  Therefore,  the  variance  of  the 
estimator  is  closely  given  by 


2 2 2 2 2 2 

1 + 2esrk  + e 1 + «?k  + 2« 

pa  — : — . — , ---  + 


4K(2Tk  + l) 


K 


9 K-l 

2Pk  E 

r=l 


(5.24) 


Therefore,  from  (5.14)  and  (5.24)  we  can  express  the  esti- 
mate of  the  error  rate  of  the  ktb  diversity  branch  acting  alone  by 

by 

Pk  “ Pk(1+t°k)  (5.25) 
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where 


<°k 


+ 


2 

e 


Substituting  (5.25)  into  the  pression  for  the  error  rate 
estimate  of  (5.5)  gives 


(5.26) 


A 

One  may  readily  compute  the  raean  and  variance  of  P 
small  co^ 


(5.26).  For 


(5.27) 


for  independent  diversity  channels. 


The  estimator  of  the  error  rate  is  biased.  For  the  systems 
of  interest,  the  effect  of  this  bias  should  be  assessed  and,  if 
the  bias  is  unacceptable,  the  filter  parameters  of  the  magnitude- 
squared  channel  transfer  function  estimator  of  Figure  4.1  can  be 
chosen  to  reduce  this  bias  for  the  expected  range  of  F^  1 s . 


From  (5.26)  and  (5.27)  the  variance  of  the  error  rate  esti- 
mate for  small  is  given  by 


(5.28) 


Examples  are  presented  in  Section  5.1.5  to  illustrate  the 
effectiveness  of  the  error  rate  estimator  analyzed  above. 
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While  the  above  calculations  of  error  bias  and  variance 
assumed  Gaussian  noise,  it  appears  that  the  results  will  not 
differ  significantly  for  atmospheric  noise.  (Calculations  of 
error  rate  with  atmospheric  rather  than  Gaussian  noise  have  been 
carried  out  by  Bello  [5.17].)  This  may  be  deduced  by  noting  that 
the  curves  of  error  rate  vs.  SNR  for  atmospheric  noise  (see 
Figures  2 and  3 of  [5.17])  are  similar  to  those  for  Gaussian  noise. 
These  comments  apply  to  all  the  error  rate  estimation  techniques 
and  will  not  be  repeated. 

5.1.3  Use  of  Diversity  Combiner  Output  Instantaneous  SNR 
Estimates  for  Error  Rate  Prediction 


In  this  section  we  address  the  problem  of  estimating  the 
error  rate  for  a flat  fading  channel  by  averaging  instantaneous 
error  rate  estimates  determined  after  pseudo-diversity  combining. 
Figure  5.3  is  a functional  block  diagram  of  the  error  rate  esti- 
mator to  be  considered.  The  estimation  procedure  is  as  follows. 
The  data  signal  is  transmitted  over  the  diversity  channel,  picked 
off  at  RF  (or  IF)  in  the  receiver,  and  sent  to  a magnitude-squared 
channel  transfer  function  estimator  of  the  type  analyzed  in  Sec- 
tion 4.1.  The  output  of  this  estimator  is  sampled  and  combined 
with  an  estimate  of  the  noise  power  to  form  an  estimate  of  the 
instantaneous  SNR.  This  instantaneous  SNR  estimate,  along  with 
instantaneous  SNR  estimates  from  the  other  diversity  branches,  is 
sent  to  a simulated  diversity  combiner.  This  simulated  combiner 
produces  an  estimate  of  the  combiner  output  SNR  at  that  instant. 

An  instantaneous  error  rate  is  then  computed  and  averaged  to  give 
an  error  rate  estimate. 

For  analysis  purposes,  the  modem  is  chosen  to  be  a two-phase 
DCPSK  system  employing  order  maximal  ratio  combining.  A 

statististically  stationary  channel  with  independent  fading  on 
diversity  branches  will  be  assumed  to  simplify ' calculation. 
However,  it  should  be  noted  that  the  method  does  not  require  any 
assumption  about  the  dependence  or  statistics  of  the  fading.  It 
requires  only  that  the  additive  nod ses  be  Gaussian  and  independent . 

For  ideal  maximal  ratio  combining  and  independently  fading 
diversity  branches,  the  instantaneous  oNR  at  the  combiner  output 
is  the  sum  of  the  input  SNR's  [5.1].  Hence,  the  SNR  at  the  com- 
biner output  and  at  the  pth  sampling  instant  will  be  estimated  by 


M 

- E 

k=l 


kp 


(5.29) 
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Bias  Corrected 


t*Vi 

where  are  the  estimates  of  the  SNR  on  the  kcn  diversity  branch 
and  at  the  pfh  sampling  instant. 

Furthermore,  recalling  [from  (5.1)]  the  expression  for  the 
error  rate  of  a two-phase  DCPSK  system  in  Gaussian  noise,  we  can 
form  an  estimate  of  what  the  ensemble  error  rate  would  be  at  the 
pfch  sampling  instant: 


J(p) 


(5.30) 


An  estimate  of  the  average  error  rate  can  be  obtained  by 
averaging  these  instantaneous  error  rates.  Hence,  the  average 
error  rate  estimate  is  given  by 
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• E 

k=l 


^kp 


(5.31) 


Representing  the  estimated  instantaneous  SNR’s  by  Eq.  (5.6), 
the  error  rate  estimate  becomes 
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a+ckP)+5kp 


(5.32) 


where  and  fi^p  are  errors  in  estimating  y^p* 

The  effectiveness  of  the  above  error  rate  estimate  is  deter- 
mined by  evaluating  its  bias  and  rate  of  convergence.  We  will 
assume  that  the  errors  are  independent  of  each  other  and  of  y^p. 
Furthermore,  the  diversity  branches  will  be  assumed  to  fade 
independently  in  the  evaluation  of  estimator  performance.  Note, 
however,  that  the  technique  itself  is  valid  for  dependent  diver- 
sity channels  since  conditional  error  rates  are  computed  only 
with  the  diversity  combined  SNR.  Therefore,  we  can  write 


■ 


_1_ 

2K 


K M 

s IT 

p=1  k-i 


-a 


lf<e 


kp 


H 


-ykp(i+tkp)i 


(5.33) 
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Using  (5.7)  to  represent  the  moments  of  the  errors,  then  in 
a manner  similar  to  that  shown  in  Section  5.1.2,  it  can  be  shown 
that  E-f f’g j-  can  be  closely  approximated  by 


{U-  1+f  £*£+*'* 


(5.34) 


where  higher  than  second-order  moments  of  the  error  were  neglected 
and  where,  for  the  system  analyzed,  the  error  rate  that  we  are 

1 

^ - .i* • _ j i *- 


estimating  is  given  by 


2 ff  <rk  + « 

k=l 


Comparing  (5.27)  with  (5.34),  we  note  that  the  biases  for 
the  two  error  rate  estimators  considered  in  this  section  are 
identical. 


The  second  moment  of  P is  found  by  squaring  and  averaging 
(5.32)  to  give  e ' “ 
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S y . (1+e.  )+y.  (1+e.  )+6.  +6. 

k=l  kP  kp'  kq  kq'  kp  ukq 


(5.35) 


Assuming  that  the  errors  in  estimating  the  instantaneous  SNR's 
at  different  sampling  instants  are  independent,  we  have 


--Si : ff  B|«"2#kpj  Eje'2ykp(1+M 

e 4lT  p-1  k'2x  [ t M ) 
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(5.36) 
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Performing  the  above  averages  and  neglecting  higher  than 
second-order  moments  of  the  error,  it  follows  that  we  can  closely 
approximate  the  second  moment  of  the  error  rate  estimate  by 


+ 


9 9 9 

1 + 2Me  + ef  E r. 

S k~l 
4K2 


(5.37) 


where 


is  defined  by 


(5.16)  and  where 


a B 

W ■ rf 

k=l. 


<rk+1)‘ 


(rk+i)2-  (Ppqrk)2 


- 1 


(5.38) 


Also,  when  M=l,  Oi(ppq)  reduces  to  $>(p  ) as  given  by  (5.20). 

In  fact,  it  is  easy  to  see  that  when  the  error  rate  esti- 

mates given  by  (5.5)  and  (5.31)  are  identical,  as  are  the  tech- 
niques given  by  Figures  5.1  and  5.3. 

For  a stationary  scatter  channel,  the  double  sum  of  (5.37) 
can  be  reduced  to  a single  sum  by  defining  p(r)  as  in  the  previous 
section.  This  gives  for  the  variance  of  the  estimate 
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where  the  terms  that  fall  off  faster  than  ~ and  higher  than 
second-order  error  moments  were  neglected. 

For  the  usual  case  of  not  strongly  correlated  diversity 
branches,  the  error  rate  computed  at  the  diversity  combiner  out- 
put will  be  very  much  lower  than  the  error  rate  of  each  branch 
acting  alone.  As  with  an  error  counting  procedure,  the  error 
rate  estimate,  variance  varies  as  the  reciprocal  of  the  error 
rate  that  it  is  desired  to  measure.  Thus,  the  variance  of  the 
Pe  estimate  for  method  2 will  normally  be  much  worse  than  for 
method  1. 

A variant  of  method  2 was  presented  in  Section  3 which  re- 
duces greatly  the  variance  of  the  Pe  estimate.  Briefly,  the 
apparent  SNR  is  reduced  by  a factor  r in  the  conditional  error 
rate  estimates.  Because  the  apparent  SNR  is  much  less,  the 
apparent  error  rate  is  higher  and  a much  smaller  variance  of  the 
error  rate  estimate  (and  bias)  is  achieved.  The  correct  error 
rate  is  estimated  by  extrapolation  with  the  aid  of  a theoretical 
error  rate  vs.  SNR  curve.  It  was  shown  in  Section  2 that  the 
procedure  is  exact  at  high  SNR  and  will  work  for  a wide  class  of 
fading  statistics,  additive  noises,  modulation  techniques,  and 
diversity  combiners. 

5.1.4  Error  Burst  Estimation 

In  this  section  we  examine  the  problem  of  estimating  error 
bursts  from  estimates  of  the  instantaneous  error  rate.  The  error 
burst  estimates  can  be  used  to  distinguish  error  bursts  caused  by 
the  channel  from  those  caused  by  the  equipment. 

The  estimator  given  by  Figure  5.3  will  be  examined  for  its 
ability  to  estimate  short-term  error  bursts.  A two-phase  DCPSK 
system  employing  maximal  ratio  combining  will  be  used  for  the 
analysis.  For  this  system,  the  error  rate  estimate  is  given  by 
[repeating  (5.32)] 

- i K "k?iykp(1+Ekp)+6kp 

Pe  " S £ e «•«>) 

P = 1 

where  M is  the  order  of  diversity,  y^.  is  the  SNR  on  the  k*"*1 
diversity  branch,  at  the  pt“  sampling  instant,  and  and  6]^ 
are  errors  in  estimating  y^  due  to  noise,  data,  ancrthe 
selectivity  of  the  channel. 
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It  will  be  assumed  that  the  samples  are  taken  close  enough 
together  such  that  the  SNR  can  be  considered  to  be  constant  in 
each  sampling  interval.  Since  the  fading  rate  of  the  channel  is 
much  slower  than  the  data  rate,  many  data  bits  will  be  received 
between  samples  of  the  output  of  the  magnitude-squared  channel 
transfer  function  estimator.  We  will  assume  that  N data  bits  are 
received  for  each  SNR  sample.  Therefore,  since,  from  (5.40), 
there  are  K SNR  samples  on  each  branch,  NK  data  bits  are 
received  in  the  time  interval  over  which  we  are  estimating  error 
bursts . 


The  fraction  of  data  bits  received  in  error  during  the  ob- 
servation interval  is  given  by 


E E 

>=i  6*1 


j3(p,£) 


(5.41) 


where 


fl(p,4) 


|1,  if  the  data  bit  in  the  pfc^ 
I sampling  interval  is  in  error 

(0,  otherwise 


In  this  section,  we  will  examine  the  use  of  the  error  rate 
estimate  given  by  (5.40)  as  a means  of  estimating  given  by 
(5.41).  The  effectiveness  of  this  estimator  is  determined  by  the 
mean-squared  difference  given  by 


MSD  = E IP  -P 


(5.42) 


By  direct  substitution,  MSD  can  be  expressed  as 
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To  evaluate  the  MSD,  we  will,  first  find  the  MSD  conditioned 
upon  ■ ykp,  6^}.  This  conditional  mean-squared  difference 


is  given  by 


, K K N N ( ) 
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Recalling  the  assumption  that  the  SNR  estimates  are  taken 
close  enough  together  such  that  the  SNR  can  be  assumed  constant 
in  a sampling  interval  and  also  recalling  (5.1),  then  for  the  two- 
phase  DCPSK  system  considered,  the  required  expectations  are  given 

by 
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Substituting  (5.45)  and  (5.46)  into  (5.44)  gives  for  the 
conditional  mean-squared  difference 
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(5.47) 


Using  the  joint  probability  density  function  for  y^p  an(*  7^ q 
given  by  (5.17),  the  conditional  mean-squared  difference  given 
above  can  be  averaged  over  the  instantaneous  signa 1-to-noise  ratio 
samples  to  give  the  mean-squared  difference  conditional  upon  the 
instantaneous  SNR  estimation  errors.  Performing  this  average 
gives 
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(5.48) 


th 

where  F^.  is  the  mean  SNR  on  the  k diversity  branch  and  p is 
the  channel  correlation  coefficient  defined  by  (5.16).  ThE%iean- 
squared  difference  can  be  closely  approximated  by  the  first-  and 
second-order  moments  of  and  6^  * Therefore,  averaging  (5.48) 
over  the  errors  in  estimating  the  Magnitude  squared  of  the  chan- 
nel transfer  function  gives  for  the  mean-squared  difference 
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where 

pe  = m 

2 rf  (rk+i> 
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is  the  average  error  rate  of  the  system. 


It  should  be  noted  that  as  K->®,  the  mean-squared  difference 
given  by  (5.49)  goes  to  zero.  However,  from  (5.34),  the  estimate 
given  by  (5.40)  is  a biased  estimate  of  error  rat£  and,  since  P^, 
given  by  (5.41)  converges  to  the  error  rate,  then  Pe  is  also  a 
biased  estimate  of  P^.  This  bias  can  be  included  in  the  expres- 
sion for  the  MSD  to  obtain 
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The  performance  of  the  error  burst  estimator  considered  in 
this  section  will  be  discussed  in  the  following  section  along 
with  the  performance  of  the  error  rate  estimators  previously 
considered . 
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In  this  section,  the  performance  of  the  estimation  techniques 
analyzed  in  the  previous  sections  will  be  presented.  For  illus- 
trative purposes,  the  channel  models  (HF,  LOS,  Tropo,  and 
Satellite)  used  in  Section  4.1  to  analyze  the  performance  of  the 
magnitude -squared  channel  transfer  function  estimator  will  also 
be  used  in  the  examples  of  this  section. 


From  (5.24),  (5.28),  and  (5.39)  the  determination  of  the 
error  rate  estimator  performances  requires  the  functional  descrip- 
tion of  the  channel  correlation  function.  To  illustrate  the 
performance  of  these  estimators,  two  representations  for  p( r) 
were  used;  one  is  the  Gaussian  correlation  function  and  is 
given  by 


PG(r)  = e 
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and  the  second  is  the  double  pole  correlation  function  given  by 


PDP(r) 


[1  + 7T  j r 1 B At]  e 
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where 


B is  the  channel  rms  Doppler  spread 
At  is  the  time  between  samples 


It  should  be  reiterated  that  the  above  representations  for 
the  channel  correlation  function  are  for  illustrative  purposes 
only,  and  are  not  meant  to  represent  the  actual  correlation 
functions  for  the  various  channels.  It  was  found  in  the  examples 
that  the  performance  of  the  error  rate  estimators  is  relatively 
invariant  for  the  above  two  correlation  functions  and,  thus,  the 
performance  of  the  estimators  may  be  robust  with  respect  to  the 
shape  of  the  correlation  function.  The  performance  presented  is 
for  the  double  pole  correlation  functio. . 

Figure  5.4  and  5.5  present  the  error  rate  estimator  perfor- 
mance for  the  estimator  analyzed  in  Section  5.1.2,  while 
Figures  5.6  and  5.7  illustrate  the  performance  of  the  estimator 
analyzed  in  Section  5.1.3.  The  estimation  errors  cr£  given  by 
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Figure  5.4  Performance  of  Error  Rate  Estimator  1 for 
HF  Channel  (10-minute  Estimation  Time) 
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(5.28)  and  (5.39)  are  used  as  the  measure  of  the  estimator  per- 
formance. The  performance  of  the  estimators  for  other  time 
intervals,  still  large  such  that  the  sums  in  (5.24)  and  (5.39) 
do  not  change  appreciably,  can  be  determined  from  these  figures 
by  noting  that  a*  decreases  as  the  inverse  of  the  square  root 

*0 

of  the  number  of  samples.  Thus,  if  the  performance  for  a time 
interval  of  four  times  the  specified  time  is  desired,  it  can  be 
determined  by  decreasing  crt  by  a factor  of  2 from  the  perfor- 

mance  presented. 

Table  5-1  presents  the  parameters  used  to  evaluate  the  error 
rate  estimation  techniques.  As  pointed  out  in  Section  2,  the  LOS 
channel  should  be  analyzed  on  a quasi-stationary  basis.  The  LOS 
channel  has  been  included  in  Table  5-1  to  make  clear  that  reliable 
average  error  rate  measurements  on  LOS  links  are  meaningless 
because  of  the  long  measurement  time  required  — 1667  minutes. 

From  Figures  5.4  to  5.7,  we  note  that  the  estimator  of  Sec- 
tion 5.1.2  is  superior  to  the  estimator  of  Section  5.1.3.  Since 
the  bias  of  Error  Rate  Estimators  1 and  2 can  be  significant  for 
HF  channels,  it  is  shown  on  these  performance  curves.  This  bias 
was  calculated  using  estimator  filter  parameters  for  which  the 
bias  Is  minimized  at  each  SNR.  Therefore,  this  represents  an 
upper  bound  on  performance.  Due  to  this  bias,  these  estimators 
are  not  recommended  for  use  on  HF  channels . For  the  other  chan- 
nels of  interest,  the  bias  can  be  made  negligible  by  proper  filter 
parameter  selection.  Use  of  the  error  amplification  technique 
presented  in  Section  3 will  also  reduce  the  error  rate  bias 
significantly . 

Figures  5.8  -5.11  present  the  performance  of  the  error  burst 
estimation  technique  analyzed  in  Section  5.1.4.  A short  time  in- 
terval (15  seconds)  is  used  in  the  example  because  good  short-term 
error  burst  estimates  can  aid  an  operator  in  quickly  distinguishing 
an  error  burst  due  to  the  channel  fading  from  one  due  to  equipment 
failure. 

From  the  figures  we  note  that  the  best  error  burst  estimates 
(small  MSD)  are  obtained  for  the  channels  for  which  good  estimates 
of  the  instantaneous  SNR  are  possible  (LOS,  Satellite,  Tropo). 

For  the  channels  that  do  not  allow  good  estimation  of  the  instan- 
taneous SNR,  the  error  burst  estimator  of  Section  5.1.4  is  of 
questionable  usefulness. 
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TABLE  5-1 

PARAMETERS  USED  TO  EVALUATE  ERROR  RATE  ESTIMATION  TECHNIQUES 


DIMENSIONLESS 


I I I I i ! I ■ I : • i ■ I .i  1 ' I 

[At  = 0.174,  ? = 0.003,  €g  = 0.0] 

DATA  RATE  - 4 Mb/s 

MSD  is  the  Mean  Squared  Difference 


'nwrn,‘ 


5 . 2 Flat  Fading  Error  Rate  Estimation  from  Measured  Fade 
Probabilities 


5.2.1  Introduction 

For  a flat  fading  Rayleigh  fading  channel,  the  probability 
that  the  signal-to-noise  ratio  fades  below  a level  can  be  related 
to  the  average  error  rate  of  the  channel.  In  this  section,  we 
use  measured  fade  probabilities  as  a means  of  estimating  the 
error  rate  of  a flat  fading  channel.  Two  techniques  for  esti- 
mating the  error  rate  are  considered.  The  first  technique  uses 
the  fade  probabilities  to  form  estimates  of  the  error  rate  each 
diversity  branch  would  have  if  acting  alone,  and  then  uses  these 
estimates  of  branch  error  rates  to  estimate  the  error  rate  of 
the  channel.  The.  second  technique  uses  the  fade  probabilities  of 
the  output  of  the  diversity  combiner  to  estimate  the  average  error 
rate.  While  simpler  to  implement  than  the  techniques  described 
in  Section  5.1,  these  techniques  are  approximate,  requiring  that 
the  SNR  be  large  compared  to  unity. 

5.2.2  Relationship  Between  Fade  Probabilities  and  Error 
Rate 


For  a nonselective  Rayleigh  fading  channel,  the  error  rate 
for  many  modems  can  be  directly  related  to  the  probability  that 
the  signal-to-noise  ratio  fades  below  a level.  In  this  section, 
we  will  determine  this  relationship. 

The  error  rate  for  many  modems  operating  over  a nonselective 
Rayleigh  fading  channel  can  be  approximated  by  [see  Ref.  5.1, 
chapter  10] : 


P « 
e 


k=l 


r.  » 1 

k 


(5.53) 


where  M is  the  order  of  diversity  and  rk  is  the  mean  SNR  on  the 
diversity  blanch.  Equation  (5.53)  is  a good  approximation 
when  rk»  1.  The  constant  S is  dependent  upon  the  order  of  diver- 
sity, the  type  of  combining,  and  the  modulation  technique. 

Table  5-2  gives  values  of  S for  some  modulation  schemes  and 
diversity  combiner  techniques. 

The  probability  that  the  SNR  at  the  output  of  the  diversity 
combiner  is  below  a level  is  approximately  given  by  [see  Ref.  5.1, 
chapter  10] 
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TABLE  5-2 


i 


[ 

i 
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VALUES  OF  S SUCH  THAT  P m 

e 


S 


IT  rk 


k=l 


Combining 

^X.J£echnique 

Modulation''1'''*. 

Scheme 

Maximal 

Ratio 

Combining 

Selection 

Combining 

Equal 

Gain 

Combining 

Coherent 

PSK 

(M-|)! 

W 

M! 

2 Jv 

M!2 

i 

Coherent 

FSK 

2M"1(M  -|>! 

2M‘1(M  --)! 

rf* 

M!  Jv 

2 M! 

Differentially 

Coherent 

PSK 

1 

2 

M! 

2 

2(M-|)! 

Noncoherent 

FSK 

2m-i 

2M‘1  M! 

vF1  Jn 

2(M-|)! 

! 

th 

rk  is  the  mean  SNR  on  the  k diversity  branch 
M is  the  order  of  diversity 
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A 


yM 

Prob . jy  s£  xj>  » 

iir k 

k=*l 


for  X « I\  , k=l , . . . ,M 

<X>1) 


(5.54) 


where 


for  maximal  ratio  combining 


= 1 for  selection  combining 


M!  (M-|>! 


■z — for  equal  gain  combining 

J-  \ I 


The  fade  probability  given  by  (5.54)  is  a good  approximation 
when  X«Tk.  Combining  Eqs.  (5.53)  and  (5.54),  it  follows  that 


Prob 


• {y*x} 


; for  X«Tk,  k=l,...,M  (5.55) 

(X>1) 


From  the  expressions  for  V and  S,  it  is  easy  to  show  that 
for  the  three  combining  techniques  considered,  the  ratio  ^ is  the 
same  for  all  combiners.  In  particular,  for  DCPSK  and  noncoherent 
FSK,  the  fade  probability  is  given  by 


Prob . jy  s x} 


2a1  X11  P 


; for  X«Tk,  k=l,...,M  (5.56) 

(X>  1) 


where  a ^ 1 for  DCPSK,  a.-^  for  noncoherent  FSK,  For  the  coherent 
detection  techniques 


Prob 


.{y  ^X} 


2 Vtt  /3M  X1^  P 


(“■I)* 


; for  X«Fk>  k=l,...,M  (5.57) 

(X>  1) 


where  /3  * 1 for  PSK,  — f°r  FSK. 
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Therefore,  the  probability  that  the  SNR  at  the  output  of  the 
combiner  falls  below  a level  is  directly  related  to  the  proba- 
bility of  error.  In  Sections  5.2.3  and  5.2.4,  we  will  examine 
two  techniques  for  estimating  error  rate  from  measured  fade 
probabilities . 

5.2.3  Use  of  Fade  Probabilities  on  Each  Diversity  Branch 
for  Error  Rate  Estimation 

In  this  section,  we  will  use  error  rate  estimates  for  each 
diversity  branch  to  predict  the  system  error  rate.  The  error 
rates  of  the  diversity  branches  are  estimated  from  their  fade 
probabilities.  Figure  5.12  is  the  block  diagram  of  the  estimator 
to  be  analyzed.  The  estimation  procedure  is  as  follows.  The 
data  signal  is  transmitted  over  the  diversity  channel,  picked  off 
at  RF  (or  IF)  in  the  receiver,  and  sent  to  the  magnitude -squared 
channel  trarisfer  function  estimator  of  the  type  analyzed  in 
Section  4.1.  The  output  of  this  estimator  is  sampled  and  com- 
bined with  an  estimate  of  the  noise  power  to  form  an  estimate  of 
the  instantaneous  SNR.  This  instantaneous  SNR  estimate  is  then 
compared  with  a threshold,  quantized  and  summed  to  form  an  esti- 
mate of  the  error  rate  on  the  kt^1  diversity  branch.  This  error 
rate  estimate,  combined  with  error  rate  estimates  from  the  other 
diversity  branches,  is  then  used  to  estimate  the  average  error 
rate. 

To  estimate  the  error  rate  for  one  diversity  branch  acting 
alone  from  fade  probabilities,  either  (5.56)  or  (5.57)  should  be 
used  with  M=l.  Therefore,  the  ensemble  error  rate  on  the  kth 
diversity  branch,  assuming  the  other  diversity  branches  are  not 
there,  can  be  approximately  related  to  the  fade  probability  by 

A Prob.jy,  £ x} 

\ ^ 1 (5.58) 

ttl 

where  y^.  is  the  instantaneous  SNR  on  the  k diversity  branch  and 
A can  be  found  from  (5.56)  or  (5.57)  for  the  modems  considered. 

In  implementing  the  technique  to  estimate  the  error  rate  on 
the  kth  diversity  branch  acting  alone,  the  fade  probabilities  will 
be  estimated  from  the  instantaneous  SNR  estimates  by 

?rob.{yk  = x}  4 E b 0.59) 

P = 1 ^ 
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Figure  5.12  Block  Diagram  of  Error  Rate  Estimator 


where 


1,  ifykp*x 


0,  otherwise 


th 


and  y^  is  an  estimate  of  y^  at  the  p sampling  instant. 

th 

Therefore,  the  error  rate  on  the  k diversity  branch  will 
be  estimated  by 


K 


& 2 


(5.60) 


k=l 


To  evaluate  the  effectiveness  of  the  error  rate  estimator 
given  in  Figure  5.12  to  estimate  error  rate,  we  must  determine 
the  convergence  of  P^.  to  P^.  From  Eq.  (5.6),  for  the  magnitude- 
squared  transfer  function  estimator  of  Section  4.1,  the  estimate 
of  the  instantaneous  SNR  can  be  expressed  as 


Xkp  “ rkp(1  + tkP)  + 6kP  (5-61) 

til 

where  y.  is  the  actual  signal-to-noise  ratio  on  the  k diversity 
branch  alid  at  the  sampling  instant.  The  parameters  e.  and 
6^p  are  errors  in  estimating  y^p  due  to  noise,  data,  and  the 
selectivity  of  the  channel.  * 

A 

From  (5.59)  and  (5.60),  the  expected  value  of  can  be 
given  by 


Prob. 


k 


x-6kP) 

1+ekp| 


(5.62) 


Using  (5.58)  to  express  the  above  probability,  averaging 
over  the  error  parameters  gives 

E{Pk}  - Pk(1  + t2)  (5.63) 

where  the  moments  of  the  errors  in  estimating  y,  are  given  by 
(5.7)  and  higher  than  second-order  terms  of  these  errors  are 
considered  negligible. 
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Therefore,  the  estimate  of  the  error  rate  on  the  k*-*1  diver- 
sity branch  is  biased  due  to  the  error  in  estimating  the 
magnitude  squared  of  the  channel  transfer  function.  The  second 
moment  of  can  be  found  from  (5.60)  to  be  given  by 


K K 


TO  “-T2  E E Ebb 

1 X K (p-lq-1  P 


(5.64) 


From  the  definition  of  b , and  noting  that  b“j  = b , it  imme- 
diately follows  that  P 


E{$ 


a o .2  K K ( X - 6, 

A Pt(l  + «2)+-4-o  E E Prob.  ytn^,y.  * 


X - 6 


kq 


KX  k 


2 2 

XT  p=lq=l 


kp'1  + <kp  kq  1 + fkq 

(5.65) 


With  the  joint  density  of  the  SNR  samples  p (y,  ,y,  ) given 

by  (5.17),  we  have  P ^ 

x-6kP  x-6ka 

l X-5k  X-6.  ) >’> 

Prob’jrkpST+7^’ykq£  1+7^1  “ J J P^-V^kp^k, 

1 r 0 0 

(5.66) 


For  analysis  purposes,  we  will  assume  that  the  SUR  samples 
are  spaced  far  enough  apart  in  time  such  that  for  all  y.  a X, 
y^  £ X,  we  have  P 


^Ppq  V^kp^kq 

Fk(l-Ppq) 


small 


small 


(5.67) 
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With  the  above  assumptions,  then  p 
approximated  by 


(ykp,ykq^  can  be  closely 


P(’'kp'ykq)  “ 


ykp’ykq  s X 


(5.68) 


Since  in  order  for  the  error  rate  to  be  related  to  the  fade 
probability  we  must  have  X«r^,  then  the  small  argument  approxi- 
mation of  (5.67)  will  be  valid  if  p is  not  too  near  1.  Substi- 
tuting (5.68)  into  (5.67),  integrating  with  respect  to  y,  and 
y,  , and  averaging  over  the  errors  give  for  the  joint  fade 
prSbability 


X“6kp 

Prob-^kPsr+7'^>y-  s 


X-  6 


) X2(l  + 2e2) 


(5.69) 


Substituting  (5.69)  into  (5.65)  gives  after  some  simplifica 

tion 


21  „ A. 

KX 


0 P2(l  + 2e2)  K K 

p,.(i + €2)  +-S — o E Z 


K p=l  q=l  1 - p 

pA  pq 


(5.70) 


For  a stationary  scatter  channel  and  equal  spacing  between 
samples,  we  can  express  p in  terms  of  r = p-q.  Therefore, 
letting  p(r)=p  in  (5.70;,  the  double  sum  can  be  replaced  by 
a single  sum  toFgive  for  the  variance  of 


2 A 2 pva+2c2)  2P2(l  + 2f2)  fcl 

yii;  pka  + f > K + ^2 r5x  *tp(r)](K-r) 

k 

(5.71) 


where 

«/)[p(r)]  = — P-  2^'  ■ 
1 - p (r) 
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If  0[p(r)]  contributes  significantly  to  the  sum  only  when 
r«K  and  if  (as  should  be  the  case)  the  second  term  is  negli- 
gible, then  the  variance  of  ^ is 


i Ik 

KX 


(1  + e ) + 


2P2(l  + 2e2)  Krl 


K 


0[p(r)] 


(5.72) 


From  (5.60)  and  (5.72),  the  estimate  of  the  error  rate  on 
the  diversity  branch  acting  alone  can  be  expressed  as 


K - V1+nk> 


(5.73) 


where 


Ek}  ■ <2 


Proceeding  as  in  (5.25)  to  (5.28),  the  mean  and  variance  for 
estimating  the  error  rate  with  the  technique  of  Figure  5.16  is 


P = P (1  + Me2) 
e ev 


cr 


2 


P 

e 


M 


E 

k=l 


(5.74) 


The  performance  of  this  estimator  will  be  illustrated  in 
Section  5.2.5.  However,  it  should  be  noted  that  in  implementing 
the  error  rate  estimator  analyzed  in  this  section,  a choice  of 
the  level  must  be  made.  From  Section  5.2.2,  the  level  on  each 
diversity  branch  must  be  much  less  than  the  mean  SNR  of  that 
branch.  This  allows  us  to  relate  the  fade  probability  on  that 
branch  to  the  error  rate  of  that  diversity  branch  acting  alone. 
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Furthermore,  from  (5.72)  the  rate  of  convergence  of  the  estimator 
Increases  as  the  level,  X,  increases. 

5.2.4  Use  of  Combiner  Output  Fade  Probabilities  for  Error 
Rate  Estimation 


5.2.4. 1 Introduction 

As  was  shown  in  Section  5.2.2,  for  high  SNR  the  probability 
that  the  instantaneous  SNR  at  the  output  of  the  diversity  combiner 
is  below  a level  is  directly  related  to  the  average  probability  of 
error.  Figure  5.13  is  a block  diagram  of  an  error  rate  estimator 
that  utilizes  the  above-mentioned  relationship.  Up  to  the  esti- 
mation of  the  instantaneous  SNR  at  the  diversity  combiner  output, 
the  estimation  procedure  is  identical  to  that  of  Error  Rate 
Estimator  2 discussed  in  Section  5.1.3.  However,  instead  of 
estimating  the  instantaneous  error  rate,  the  error  rate  estimator 
of  Figure  5.17  compares  an  estimate  of  the  instantaneous  SNR  at 
the  diversity  combiner  output  with  a level  to  form  an  estimate  of 
a fade  probability.  These  fade  probability  estimates  are  then 
used  to  estimate  the  average  probability  of  error  of  the  channel. 

From  Eqs.  (5.56)  and  (5.57),  the  error  rate  of  the  systems 
considered  can  be  estimated  by 


B Prob . jy  <:  x| 

7 " 


(5.75) 


where  B is  a constant  dependent  upon  the  particular  modulation 
scheme  used,  M is  the  order  of  diversity,  and  y is  the  instanta- 
neous SNR  at  the  output  of  the  diversity  combiner.  The  proba- 
bility that  the  SNR  at  the  combiner  output  fades  below  a level  X 
will  be  estimated  by 


where 


Prob . £ X j 
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K 
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P=i 
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i>  if  yp  £ x 

0,  otherwise 


(5.76) 
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Figure  5.13  Block  Diagram  of  Error  Rate  Estimator 


A til 

and  y is  the  estimate  of  the  SNR  in  the  p sampling  interval. 

Substituting  Eq.  (5.73)  into  (5.72)  gives 


A 

P 

e 


B 


K 

E 

-Eri 


3C“ 


K 


(5.77) 


5. 2. 4. 2 Instantaneous  Combiner  Output  SNR 

A 

To  evaluate  the  effectiveness  of  Pg  given  by  (5.77)  to  esti- 
mate the  average  error  rate,  it  is  necessary  to  relate  the  esti- 
mate of  the  instantaneous  SNR  at  the  output  of  the  combiner  to 
its  actual  value.  We  will  derive  this  relationship  for  two  of 
the  combining  techniques  considered  in  Section  5.2.2.  It  will  be 
assumed  that  [see  Eqs.  (4.250),  (5.6)] 

*kp  = V1+fkp>  + V (5-78) 

til 

the  k diversity  branch  and  at  the 
is  the  estimate  of  y^  . 

For  selection  combining,  the  combiner  output  SNR  is  the 
largest  of  the  input  SNR.  With  small  errors  in  estimating  the 
magnicude  squared  of  the  channel  transfer  function,  the  diversity 
branch  with  the  largest  estimated  SNR  will  usually  have  the 
largest  instantaneous  SNR.  Therefore,  we  can  relate  the  estimated 
and  actual  SNR  at  the  output  of  the  diversity  combiner  at  the 
pttl  sampling  instant  by 


where  y.  is  the  actual  SNR  on 
pfck  sampling  instant,  and  y^p 


(1  + e ) + 6 

P P 


where 


= 0 
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(5.79) 


i-trnim  n 


iffHiifrr'f-f'Y'n'  nlinn1 


and  € and  Cg  are  the  errors  in  estimating  the  magnitude  squared 
of  the  channel  transfer  function  and  IZ  is  the  mean  SNR  of  the 
diversity  branch  with  the  largest  SNR  at  pfch  sampling  instant , 

For  maximal  ratio  combining,  the  instantaneous  SNR  at  the 
output  of  the  diversity  combiner  is  the  sum  of  the  instantaneous 
SNR's  at  the  input  to  the  combiner.  Therefore,  the  estimate  of 
the  SNR  at  the  output  of  a maximal  ratio  combiner  at  the  pfc^ 
sampling  instant  is  given  by 


M 


where  M is  the  order  of  diversity. 

Using  (5.78)  to  express  the  estimates  of  the  instantaneous 
SNR,  it  follows  that 


?P " k?i  V1+V  + 5P 


(5.81) 


where 


M 


’ Z « 


k=l 


kp 


(5.82) 


5. 2. 4. 3 Error  Rate  Estimation 


The  performance  of  the  error  rate  estimator  presented  in 
Figure  5.17  will  be  determined  in  this  section.  Since  the  fade 
probabilities  for  equal  gain,  maximal  ratio,  and  selection  com- 
bining are  approximately  related  by  (5.54),  the  error  rate  for  these 
combining  techniques  can  be  estimated  by  simulating  a selection 
combiner  and  then  adjusting  for  the  actual  combiner  employed. 

The  error  rate  for  selection  combining  and  two-phase  DCPSK 
modulation  is  found  from  (5.56)  to  be  given  by 


M!  Prob.jy  s xj 
2 XM 


(5.83) 


where  y is  the  instantaneous  SNR  at  the  output  of  the  selection 
combiner.  The  fade  probability  will  be  estimated  by 
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where 

!1  , if  ypsX 

(5.85) 

0 , otherwise 

In  (5.85),  y is  the  estimate  of  the  selection  combiner  out- 
put at  the  ptb  sampling  instant,  and  can  be  represented  by  (5.79). 
To  determine  the  performance  of  Error  Rate  Estimator  4,  we  will 
consider  the  case  of  equal  mean  SNR  on  each  diversity  branch.  For 
this  case,  the  estimate  of  the  instantaneous  SNR  at  the  combiner 
output  is  given  by 


y - y (1  + c ) + 6 (5.86) 

P P P P 

where 

EW = EW = 0 

Ebp}  - <2 
*{<#  - & 


and  where  F is  the  mean  SNR  on  each  of  the  diversity  branches. 
From  (5.83)  and  (5.84),  the  error  rate  for  selection  combining  is 
estimated  by 
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mi  y 
EL 


p 


2K 


X* 


(5.87) 


The  mean  value  of  the  error  rate  estimate  is  given  by 
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K ( X - 6 

M!  V Prob.  y p<Y+^ 


(5.88) 


Using  (5.54)  to  express  the  fade  probability,  it  follows  that 


={&s>  = 


K (/X  - 6 f 

. ” s iM 


(5.89) 


From  Table  5-2,  we  can  approximate  (5.89)  by 


rA  , , K (/X-  6 

ik 


(5.90) 


and  neglecting  higher  than  second-order  error  moments,  we  can 
write 


(5.91) 


Therefore,  the  error  rate  estimate  is  biased  due  to  the 
errors  in  estimating  the  magnitude  squared  of  the  channel  trans- 
fer function.  The  second  moment  of  the  selection  combiner  error 
rate  estimate  is  found  from  (5.87)  to  be  given  by 


K K 


<M!>  2 2 E{Vq} 

p=l  q-1  * M 


4K2  X2M 


(5.92) 
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Since  ^(bp)  *=  E|b  j , then 
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p^q 


To  evaluate  the  joint  probability  that  two  samples  of  the 
combiner  output  SNR  are.  below  some  levels,  we  first  express  this 
joint  fade  probability  as 


U.  v 

Prob. jy  £U,y  <;  v}  = f r p(y  ,y  ) dy  dv 
l/p  'q  I J0  £ ^VV  q rp  ' q 


(5.94) 


Using  the  small  argument  approximation  for  the  joint  density 
of  PCyjcpjy^q)  given  by  (5.68),  we  can  obtain  a small  argument 
approximation  for  p(y  ,yq)  by  recalling  that  for  selection 
combining  the  SNR  at  the  output  of  the  diversity  combiner  is 
equal  to  the  largest  of  the  input  SNR's.  Therefore,  it  follows 

that  M2  y^"1  y**"1 

p(yp>yq)  w (5.95) 


m- 

k=l  ' 


Thus,  we  can  write 


( X " 6 x - 6 ) 

Prob  .<  y <:  r—r — ^ , y <;  , 

) 'p  1 4-  f ’ 'a  1 4-  c 

p q) 


Mm 


(5.96) 
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Performing  the  above  expectations,  neglecting  their  second-order 
error  terms,  and  substituting  the  end  result  into  (5.93),  the 
variance  of  the  error  rate  estimate  is  approximately  given  by 


2 M!  _ 
or  *•  w Pc 

P„  2K  X**  S 


\,m±n  2, 


2p^  r 

+ 1 + M(M  + 1)  €2  + 

K L 


2 X 


M(M-  Dtgr2  1 


K-l 

E ♦j.tpCr)] 


r=l 


(5.97) 


where 


ii)1[p(r)]  = 1 


r 2 1M 

1 - P (r) 


- 1 


A stationary  scatter  channel  with  p(r) = p (where  r = ]p-q| ) was 
assumed,  where  p is  the  channel  correlation  coefficient  defined 
by  ('■.16).  Pq 

To  estimate  the  error  rate  when  maximal  ratio  diversity  com- 
bining is  used,  it  is  only  necessary  to  use  Table  5-1  to  determine 
that  the  error  rate  for  maximal  ratio  combining  is  related  to  the 
error  rate  for  selection  combining  by 


-L  P 

Ml  S 


(5.98) 


Therefore,  the  error  rate  for  maximal  ratio  combining 
(similarly  for  equal  gain  combining)  can  be  estimated  by  simu- 
lating a selection  combiner.  From  (5.91),  (5.97),  and  (5.98), 
it  immediately  follows  that 
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(5.99) 


A comparison  of  the  flat  fading  error  rate  estimator 
analyzed  in  this  section  with  the  error  rate  estimators  previously 
analyzed  is  most  conveniently  performed  by  considering  examples. 

In  the  following  section  we  will  present  examples  to  illustrate 
the  performance  of  the  estimators  analyzed  in  Sections  5.2.3  and 
5.2.4. 


5.2.5  Performance  of  Error  Rate  Estimators 

In  this  section  we  discuss  the  performance  of  the  two  error 
rate  estimators  analyzed  in  Sections  5.2.3  and  5.2.4.  For  the 
channel  models  considered  in  Section  4.1,  examples  will  be  pre- 
sented to  illustrate  the  effectiveness  of  using  fade  probabilities 
to  estimate  the  error  rate  for  a flat  fading  channel.  Since  a 
channel  correlation  function  is  required  to  evaluate  the  variance 
of  the  error  rate  estimates,  we  will  use  the  correlation  functions 
given  by  (5.51)  and  (5.52). 

The  comments  of  Section  5.1.5  regarding  the  channel  correla- 
tion function  and  the  time  interval  over  which  the  error  rate  is 
measured  also  apply  in  this  section.  However,  for  sane  of  the 
channel  models  considered,  samples  of  the  SNR  spaced  every  T sec- 
onds, where  T is  the  duration  of  the  averaging  filter  of  the 
magnitude-squared  channel  transfer  function  estimator  will  result 
in  the  small  argument  approximation  [see  Eq.  (5.67)]  not  being 
valid.  For  these  channels  SNR  samples  spaced  farther  apart  will 
be  used. 

For  the  examples,  the  level  X was  chosen  to  be  equal  to 
one-fourth  of  the  mean  SNR.  With  the  level  set  as  above,  the 
assumption  used  in  the  derivations  that  X«F  should  not  result 
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in  a large  degradation  of  the  error  rate  estimators.  We  should 
note  that  the  level  is  an  estimator  parameter  that  affects  the 
performance  of  the  estimator  and  to  achieve  the  performance  pre- 
dicted would  necessitate  prior  estimation  of  the  mean  SNR. 

Figures  5.14  and  5.15  demonstrate  the  performance  of  the 
Error  Rate  Estimators  given  by  Figures  5.12  and  5.13,  respectively. 
The  rms  Doppler  spreads  and  the  estimation  times  for  the  various 
channels  are  presented  in  Table  5-1  of  Section  5.1.  As  previously 
mentioned  in  Section  5.1,  the  performance  of  the  LOS  channel  is 
presented  for  comparison  purposes  only.  The  long  estimation 
times  required  for  the  LOS  channel  make  average  error  rate  esti- 
mation impractical. 

Comparing  Figure  5.14  with  5.15,  it  is  seen  that  Error  Rate 
Estimator  3 is  superior  to  Error  Rate  Estimator  4. 

Comparison  of  the  estimators  analyzed  in  Section  5.2  with 
those  of  5.1  will  be  delayed  until  the  last  flat  fading  error 
rate  estimator  is  presented  and  analyzed  in  Section  5.3. 

5 . 3 Flat  Fading  Error  Rate  Estimation  from  SNR  Estimates 

5.3.1  Introduction 

For  many  systems  it  is  possible  to  express  the  error  rate 
of  a flat  fading  channel  in  terms  of  the  mean  SNR  on  each  of  the 
diver - .ty  branches.  In  this  section  we  examine  the  use  of  in- 
stantaneous SNR  estimates  to  form  estimates  of  the  mean  SNR.  The 
error  rate  will  be  estimated  from  these  mean  SNR  estimates. 

Also,  by  using  the  estimate  of  the  branch  envelope  correla- 
tion coefficient  estimated  in  Section  4.3,  the  error  rate  for 
dependently  fading  diversity  branches  can  be  estimated. 

5.3.2  Flat  Fading  Error  Rate  Estimation  for  Independently 
Fading  Diversity  Branches 

The  error  rate  for  many  digital  systems  can  be  expressed  in 
terms  of  the  mean  SNR  on  each  of  the  diversity  branches . There- 
fore, for  these  systems,  estimates  of  the  mean  SNR  on  each  diver- 
sity branch  can  be  used  to  estimate  the  error  rate.  Figure  5.16 
is  the  block  diagram  of  the  error  rate  estimator  that  will  be  con- 
sidered in  this  section.  The  estimation  procedure  is  as  follows. 
The  data  signal  is  transmitted  over  the  diversity  channel,  picked 
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off  at  RF  (or  IF)  in  the  receiver,  and  sent  to  the  magnitude- 
squared  channel  transfer  function  estimator  analyzed  in  Section 
4.1.  The  output  of  this  estimator  is  sampled  and  combined 
with  an  estimate  of  the  noise  power  to  form  an  estimate  of  the 
instantaneous  SNR.  These  instantaneous  SNR  estimates  are  then 
averaged  to  form  an  estimate  of  the  mean  SNR  on  each  of  the 
diversity  branches.  This  mean  SNR  estimate,  along  with  the  mean 
SNR  estimates  on  the  other  diversity  branches,  is  used  to  form 
an  estimate  of  the  error  rate. 

Section  4.3.4  addressed  the  problem  of  estimating  the  mean 
SNR  from  the  received  data  signal.  From  (4.253)  and  (4.256),  an 
estimate  of  the  mean  SNR  on  the  diversity  branch  can  be 
expressed  as 


rk<1+V 


(5.100) 


where 


\ - 0 

9 K'1  9 

1+*p+2X>  <»> 

In  (5.100),  e [see  Eq.  (4.104)]  is  the  rms  fractional  error  in 
estimating  the  magnitude  squared  of  the  channel  transfer  function, 
K is  the  number  of  instantaneous  SNR  samples,  and  p(*)  [see  Eq. 
(5.16)]  is  the  normalized  envelope  correlation  function  evaluated 
at  the  spacing  between  samples . 

For  analysis  purposes,  we  will  evaluate  the  effectiveness 
of  this  error  rate  estimator  for  systems  whose  error  rate  can  be 
closely  represented  by 

Pe  - -a-  (5.10D 

7K 

k=l 

where  M is  the  order  of  diversity  and  S is  a constant  dependent 
upon  the  type  of  combiner  used  and  the  particular  modulation 
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scheme  employed.  Using  (5.101)  provides  a conservative  estimate 
of  estimator  accuracy  at  low  SNR.  Values  of  S for  several  common 
combining  techniques  and  modulation  schemes  are  given  in  Table  5-2 
of  Section  5.2. 

From  (5.101),  the  error  rate  can  be  estimated  by 


(5.102) 


Substituting  the  expression  for  the  estimate  of  the  mean  SNR  on 
the  diversity  branch  as  given  by  (5.100)  and  averaging  gives 


(5.103) 


where  $ was  assumed  to  be  uncorrelated  with  8 if  n^t,  terms 
that  decrease  faster  than  1/K  where  neglected,  and  only  second- 
order  error  terms  were  considered  to  be  significant. 

From  the  above  expression  for  e/p  j>,  we  note  that  this 
error  rate  estimate  is  asymptotically  unbiased.  That  is,  as 

K-”>  E{PehV 

The  variance  of  the  estimator  can  be  found  to  be  given  by 


(5.104) 


An  rms  fractional  error  for  estimating  error  rate  can  be 
given  by 


(5.105) 
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Therefore,  the  performance  of  the  estimator  does  not  depend 
upon  the  error  rate  to  be  estimated.  It  does  depend  upon  the 
number  of  SNR  samples  K,  the  order  of  diversity  M,  rms  fractional 
error  in  estimating  the  magnitude  squared  of  the  channel  transfer 
function  e , and  the  correlation  between  instantaneous  SNR 
samples . p 

Figure  5.17  illustrates  the  effectiveness  of  Error  Rate 
Estimator  5 to  estimate  the  flat  fading  error  rate.  As  mentioned 
in  Section  5.1,  the  performance  of  the  LOS  channel  is  for  compari- 
son purposes  only.  The  long  estimation  time  for  the  LOS  channel 
makes  average  error  rate  estimation  impractical. 

5.3.3  Error  Rate  Estimation  for  Correlated  Fading  on 
Diversity  Branches 

Previously,  we  have  assumed  that  the  diversity  branches  have 
faded  independently.  For  many  systems,  the  branches  do  not  fade 
independently  and  this  dependent  fading  results  in  an  increase  in 
the  error  rate  since  the  full  diversity  effect  is  not  realized. 

To  determine  the  effect  of  dependent  fading  upon  the  error 
rate,  we  note  that  Pierce  [5.3]  has  pointed  out  that  the  density 
function  of  the  instantaneous  SNR  at  the  output  of  a maximal-ratio 
diversity  combiner  can  be  given  for  small  y by 

P(y)  ” (M  -iVdcT  ^'M'1  (5'106> 


where  det  is  the  determinant  of  the  moment  matrix  of  the  set 
|t,  (f , t ) 1 of  complex  Gaussian  random  variables  and  T,  (f,t)  is  the 
time  varying  transfer  function  on  the  kt  diversity  branch. 

At  high  mean  SNR,  most  of  the  errors  are  caused  by  the 
occasional  deep  fades  of  the  signal  and  a close  approximation 
to  the  error  rate  is  obtained  by  using  (5.105).  Bello  [5.4] 
notes  that  this  procedure  yields  asymptotically  correct  results 
for  high  SNR's  and  upper  bounds  for  all  SNR's.  Using  this  approx- 
imation, we  find  that  the  error  rate  depends  upon  the  diversity 
channel  correlations  only  through  the  determinant.  Thus,  one  may 
write 

error  rate  with  correlated  fading 
error  rate  with  uncorrelated  fading 

moment  matrix  determinant  with  uncorrelated  fading 
moment  matrix  determinant  with  correlated  fading 

(5.107) 
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Figure  5.17  Performance  of  Error  Rate  Estimator  5 for 
Tropo,  Satellite,  LOS,  and  HF  Channels 


With  P„  and  Pr>  respectively,  the  error  rates  with  uncorre- 
lated and  correlated  fading  on  diversity  branches  and  with  My 
and  Mr,  respectively,  the  montent  matrix  determinants  with 
uncorrelated  and  correlated  fading,  we  have 

pc“^pu  (5’108: 

The  above  moment  matrix  determinants  are  given  by 

**0  ■ clrk 

k“l 


Mc“c  ffrk  det 
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1 ^-l.M 
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c|jT  rk|det  ^ 


where  C is  a constant  and  G^j  is  the  complex  branch  correlation 
coefficient  defined  by  (4.25$);  that  is, 


G.  . = 
kj 


E{T*(f,t)Tk(f,t)} 


E J|Tj(f  ,t)|2  E |Tk(f,t)|2 
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For  the  techniques  involving  the  use  of  the  magnitude- 
squared  channel  transfer  function  estimates,  we  are  only  able 
to  estimate  | G^j  | . The  estimation  of  pjcj  = I G^j  | 2 was  presented 
and  analyzed  in  Section  4.3.5.  As  will  be  demonstrated  below, 
the  error  rate  for  a system  with  correlated  fading  only  on 
"adjacent  branches"  is  a function  of  pk|  and,  thus,  can  be  esti- 
mated from  the  estimates  of  the  magnitude  squared  of  the  channel 
transfer  function. 

For  the  case  of  correlated  fading  only  on  adjacent  diversity 
branches,  it  follows  that 

f M ) 

Mc  - C /T  det  ^ (5.111) 

(k“l) 


where  A^  is 


an  MxM  matrix  given  by 
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(5.112) 


The  determinant  of  can  be  evaluated  by  expanding  in  terms 
of  cofactors.  This  e mansion  yields  a recursive  formula  for 
evaluating  the  determinant  of  A*..  In  particular,  it  follows  that 


det  ^ - det  Am_1 


det  *M-2 


(5.113) 


where 

det  A^  = 1 

det  A2  = 1 - | G^2 I 2 

are  used  to  initiate  the  recursive  relationship. 
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Since  det  must  be  greater  than  zero,  then  the  representa- 
tion given  by  (5.112)  implies  that  for  M large  that  |G^js  1/2. 

Of  course,  the  above  form  should  not  be  assumed  unlessJ  | (h,c  | sO.l, 
or  so.  J 

For  the  cases  of  highly  correlated  fading,  |G^|«0,9,  or  so, 
then  the  above  representation  is  not  valid.  However,  it  may  be 
possible  to  upper  and  lower  bound  the  moment  matrix  determinant. 

Figure  5.18  presents  a block  diagram  for  estimating  the 
error  rate  over  a diversity  system  with  correlated  fading  on  the 
diversity  branches. 

The  estimation  procedure  is  as  follows.  The  data  signal  is 
transmitted  over  the  diversity  channels  and  picked  off  at  RF  (or 
IF)  in  the  receivers  and  sent  to  magnitude-squared  channel  trans- 
fer function  estimators.  The  output  of  these  estimators  is  sent 
to  the  error  rate  estimator  analyzed  in  Section  5.3.2  to  obtain 
an  estimate  of  the  error  rate  assuming  uncorrelated  fading  on 
the  diversity  branches. 

The  estimates  of  the  magnitude  squared  of  the  channel  trans- 
fer function  are  also  sent  to  branch  envelope  correlation  coeffi- 
cient estimators  of  the  type  proposed  in  Section  4.3.5.  These 
branch  envelope  correlation  coefficient  estimates  are  used  to 
estimate  the  determinant  of  the  normalized  moment  matrix.  This 
determinant  estimate  is  then  combined  with  the  estimate  of  the 
error  rate  that  assumed  uncorrelated  fading  to  estimate  the  error 
rate  with  correlated  fading. 

To  illustrate  the  effectiveness  of  the  error  rate  estimator 
given  by  Figure  5.18  to  estimate  the  error  rate  for  correlated 
fading  on  diversity  branches,  an  example  will  be  presented.  For 
analysis  purposes,  a dual  diversity  system  will  be  considered. 

For  this  system  the  error  rate  with  correlated  fading  is  closely 
approximated  by 


(5.114) 


where  the  estimation  of  p*  « was  considered  in  Section  4.3.5  and 
Appendix  B.  Therefore,  tne  error  rate  with  correlated  fading 
will  be  estimated  by 


A 


(5.115) 


5-62 


Assuming  that  the  correlated  fading  on  the  diversity  branches 
does  not  significantly  affect  the  performance  of  the  error  rate 
estimator  of  Section  5.3.2,  then 


A 


Pyd+a) 


(5.116) 


where 


E{a}  ~>  0 as  K (the  number  of  samples) 


4 

K 


K-l  , 

+ S,  p ( m) 

m=l 


CO 


Fur thetmore , from  (4.267)  and  (4.268)  the  estimate  of  the 
branch  envelope  correlation  coefficient  can  be  given  by 


P 


12 


(5.117) 


where 


E{/?}  -»  0 as  J (number  of  independent  samples)  ->  00 

EK}  ‘ j(3  + 6°12  + p12  + 6pL) 


The  above  expressions  are  applicable  for  the  channels  that 
allow  an  accurate  estimation  of  |T(f,t)|2. 

Substituting  (5.116)  and  (5.117)  into  (5.115)  gives 


r (i  + «) 

p = 

C 1 - P12-  0 


(5.118) 


With  the  number  of  samples  large  enough  so  that  good  esti- 
mates of  Py  and  p| 2 can  be  obtained,  then  the  estimate  of  the 
error  rate  with  correlated  fading  can  be  accurately  expressed  in 
terms  of  the  first  and  second  moments  of  a and  j3.  Therefore,  we 
can  write 
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(5.119) 


A 


l + a 


and  it  follows  that 


3 + 6p12+p12  + 6 


(5.120) 


Thus,  the  estimate  of  the  error  rate  with  correlated  fading  on 
the^diversity  branches  is  asymptotical  ly  unbiased.  The  variance 
of  P can  be  found  from  (5.119)  to  be  given  by 


(3  + 6p12+p12  + 6pl2)! 


(5.121) 


where  a and  /?  were  assumed  to  be  uncorrelated  and  higher  than 
second-order  moments  of  a and  /S  were  considered  to  be  negligible. 


Figure  5.19  shows  the  effectiveness  of  estimating  the  error 
rate  with  correlated  fading  on  the  diversity  branches.  In  this 
figure,  it  was  assumed  that  one  independent  sample  of  |T(f,t)|2 
can  be  taken  every  1/B  seconds.  Furthermore,  for  this  figure  we 
let  K = J,  p(m)  = 0 if  m/0.  Therefore,  from  (5.121)  we  can 
approximate  the  rms  fractional  error  by 


(3  + 6p12+p12  + 6p12) 

pi2r 

BT 


1/2 


(5.122) 


In  Figure  5.19,  20  log 


BT  was  plotted. 
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ESTIMATION  SNR  = 20  log 


Figure  5 >19  Performance  of  Error  Rate  Estimator  with  Correlated  Fading 
on  Diversity  Branches 


5.3.4  Comparison  of  Flat  Fading  Error  Rate  Estimation 
Techniques 


In  Sections  5.1  to  5.3,  we  proposed  five  techniques  for  esti- 
mating the  flat  fading  error  rate  from  the  received  data  signal. 
The  bias  and  convergence  rate  for  these  estimators  were  found  and 
examples  illustrating  performance  were  presented. 

Error  Rate  Estimators  1 and  2 used  instantaneous  error  rate 
estimates  to  estimate  the  flat  fading  error  rate.  These  estima- 
tors were  found  to  be  very  sensitive  to  errors  due  to  the  selecti- 
vity of  the  channel,  which  does  not  allow  accurate  estimates  of 
small  instantaneous  SNR's.  For  these  error  rate  estimators,  the 
bias  in  estimating  error  rate  is  highly  dependent  upon  the  error- 
in  estimating  the  instantaneous  SNR  due  to  the  channel  selectivity 
tg  [see  Eqs.  (5.14)  and  (5.34)].  In  Section  3,  an  error  amplifi- 
cation technique  using  instantaneous  SNR  estimates  is  presented. 
This  technique  is  less  sensitive  to  selectivity-induced  estimation 
errors . 

Error  Rate  Estimators  3 and  4 used  fade  probabilities  to 
estimate  error  rates.  One  problem  that  ax'ose  in  implementing 
these  techniques  was  the  choice  of  a level  which  affected  the  per- 
formance of  the  estimators.  These  estimators  were  not  as  sensi- 
tive to  channel  selectivity  (this  may,  in  part,  be  due  to  the 
analysis  that  used  a fade  distribution  linear  with  the  level)  and 
converged  more  rapidly  at  high  SNR's  than  the  corresponding  tech- 
nique of  Error  Rate  Estimators  1 and  2. 

Error  Rate  Estimator  5 used  SNR  estimates  to  estimate  the 
flat  fading  error  rate.  This  technique  converged  the  fastest  of 
all  the  techniques  considered  and  was  the  least  sensitive  to 
channel  selectivity  induced  errors.  In  fact,  as  shown  in  (5.105), 
the  filter  parameters  of  the  magnitude-squared  channel  transfer 
function  estimator  should  be  chosen  to  minimize  the  rms  fractional 
error  given  by  (4.104). 

In  conclusion,  Error  Rate  Estimator  5 is  the  recommended 
technique  for  estimating  the  flat  fading  error  rate  with  indepen- 
dently fading  diversity  branches,  and  Error  Rate  Estimator  2 can 
be  used  to  estimate  error  bursts  for  channels  that  allow  accurate 
instantaneous  SNR  estimation. 
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5 . 4 Irreducible  Error  Rate  Estimation  Due  to  Time  and 

Frequency  Selective  Fading 

In  this  section,  we  address  the  problem  of  estimating  the 
irreducible  error  rate  di:>*  ;o  time  and  frequency  selective  fading. 
The  irreducible  error  rate  s the  error  rate  that  cannot  be  re- 
duced by  increasing  the  ' gnal-to-noise  ratio. 

Two  types  of  systems  will  be  considered.  The  first  irredu- 
cible error  rate  to  be  estimated  is  that  due  to  time-selective 
fading  in  differential  phase-modulated  systems.  The  second  irre- 
ducible error  rate  is  that  due  to  frequency-selective  fading  in 
FDM-FM  systems.  The  frequency-selective  fading  causes  the  appear- 
ance of  an  intermodulation  distortion  noise  at  the  discriminator 
output  which  results  in  an  irreducible  error  rate. 


5.4.1  Irreducible  Error  Rate  of  Differential  Phase 
Modulated  Systems 

Let  the  reference  vector  for  the  i*-*1  tone  be  given  by 


+ n1 


ci 


+ jn 


si 


(5.123) 


where  r|  is  the  magnitude  of  the  fading  variable  and  is  its 
angle.  The  in-phase  (cosine)  additive  term  is  given  by  n^  and 
ng-£  is  the  quadrature  (sin)  term.  Using  capital  letters  to 
represent  complex  variables,  we  can  write  for  the  reference 
vector 


V1  - R*  + N'  (5.124) 

and  for  the  received  vector 

V = R + N (5.125) 


where  the  subscripts  i are  dropped  off  to  simplify  the  notation. 
The  complex  fading  variable  R includes  the  modulated  phase  angle 
which  can  be  assumed  to  be  zero  for  the  analysis  given  below. 

For  the  fast  fading  case,  we  have 

2 

E[RR*]  = E[R'  R'*J  = S (signal  power) 

and 

E[R'  R*]  = E[RR'*]  = p e^  S2 
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(5.126) 


where  p and  B are  the  magnitude  and  phase  of  the  channel  corre- 
lation coefficient  between  successive  frames. 

For  a complex  Gaussian  process , the  effects  of  the  time 
variations  in  R can  be  modeled  exactly  by  assuming 

R « e“^2  Rq  + Z 

and 

R'  « ej^/2  Rq  + Z'  (5.127) 

where  Rn  is  the  completely  correlated  part  of  R and  R'  and  Z and 
Z1  represent  the  zero-mean  independent  parts  of  R and  R' . Thus, 
we  have 


E[RR'*]  = e“j/?  E [R0Rg]  = pe'^S2  (5.128) 

The  variance  of  the  independent  portions  of  R are  given  by 
E [Z  Z*]  = E [(R  - R0)  (R*  - Rg)] 

= S2  - 2p  S2  + p S2  (5.129) 


or 


E[Z  Z*]  = S2(l  - p) 


(5.130) 


Using  the  above  results,  we  can  reduce  the  fast  fading 
results  to  an  equivalent  slow  fading  result  by  noting  that 

V=R+N=Rq+Z+N  (5.131) 


Therefore,  we  can  define  an  equivalent  SNR  by 


r 

eq 


R.-.>  Rg 


5 _ 


E [Z  Z*]  + E[N  N*] 


(5.132) 


or 
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(5.133) 


ra 


S2(l  -p)  + 


T 


where  T is  the  pulse  duration. 

The  slow  fading  SNR  is 

r - (5.134) 

Thus,  the  equivalent  SNR  can  be  expressed  in  terms  of  the  slow 
fading  SNR  by 


, or 

eq  (1  - p)r  + 1 


(5.135) 


Since  p is  the  correlation  coefficient  between  successive 

12  2 2 

frames,  then  p is  very  near  to  1,  and  1 -p  B T . Therefore, 

the  equivalent  SNR  can  be  approximated  by 


r(i.W)+i 


(5.136) 


The  irreducible  SNR  is  the  limit  of  (5.136)  as  That  is, 


irr 


12  2 2 


(5.137) 


The  irreducible  error  rate  is  then  found  by  using  I\rr  as 
the  mean  SNR.  Figure  5.20  is  a block  diagram  of  the  irreducible 
error  rate  estimator.  The  estimation  procedure  is  as  follows. 
The  data  signal  is  transmitted  over  the  channel,  picked  off  at 
RF  (or  IF)  in  the  receiver,  and  sent  to  one  of  the  rms  Doppler 
spread  estimators  analyzed  in  Section  4.3.2.  This  estimate  of 
the  rms  Doppler  spread  is  used  to  estimate  the  irreducible  SNR, 
which  is  then  combined  with  irreducible  SNR  estimates  of  the 
other  diversity  branches  to  form  an  estimate  of  the  irreducible 
error  rate. 
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Transmitted 
Data  Signal 


Figure  5. .20  Estimation  of  Irreducible  Error  Rate  Due  to  Time  Selective  Fading 


As  an  example  of  the  irreducible  error  rate  estimator  per- 
formance, we  will  determine  the  ability  of  the  estimator  given 
by  Figure  5.20  to  estimate  the  Irreducible  error  i*ate  for  the 
case  where  the  irreducible  error  rate  is  related  to  the  irredu- 
cible SNR's  on  the  diversity  branches  by 


irr 


M 


<k> 


ff  ‘irr 
k=>l 


Therefore,  we  will  estimate  by 

vM  M 


hrr-(§^T  a < 

k~l 


(5.138) 


(5.139) 


2 2 th 

where  B^  is  the  estimate  of  on  the  kL  diversity  branch. 

We  will  consider  two  subcases  of  this  example.  The  first 
will  assume  that  all  the  diversity  branches  have  the  same  rms 
Doppler  spread  and  that  one  Doppler  spread  estimation  is  per- 
formed. For  this  case,  we  have 


(5.140) 


For  the  rms  Doppler  spread  estimators  analyzed  in  Section 
4.3.2,  we  can  express  b2  by 


(5.141) 


2 

where  E(a]  * 0 and  Eja  j depends  upon  the  rms  Doppler  spread  esti- 
mations technique  employed.  Thus,  the  irreducible  error  rate 
estimate  can  be  represented  by 


(5.142) 


Assuming  that  a sufficiently  long  enough  estimation  time  is 
used  such  that  the  irreducible  error  rate  estimate  can  be 
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accurately  approximated  by  the  first-  and  second-order  moments 
of  a,  then 


l + M(M  - 1)  h“2} 
2 ^ 

e^2}J 

(5.143) 


i 2i  2 

Since  Eja  f decreases  as  the  time  used  to  estimate  B in- 
creases, then  Eia^T  will  probably  be  much  less  than  1 

it  follows  that L 


(5.144) 


For  the  rms  Doppler  spread  estimation  techniques  considered 
in  Section  4.3.2,  the  estimate  of  is  biased.  Therefore,  the 
irreducible  error  rate  estimate  will  also  be  biased. 

The  variance  of  the  irreducible  error  rate  estimate  is  given 
by 


or. 


- S' 


irr 


/i  9 9\2M  or  /'*>-]2M-2  „ 

T J HZ|e{b^}J  e{«  } (5.145) 


Recalling  that  E-jB*“|  and  Eja  } are  the  mean  and  variance  for 
estimating  B , then  it  should  be  noted  that  the  effectiveness  of 
estimating  P^rr  depends  upon  the  rms  Doppler  spread  estimation 
technique  used. 


The  second  example  to  be  considered  using  (5.139)  is  the 
case  in  which  the  rms  Doppler  spread  on  each  branch  is  estimated 
separately.  For  this  example  we  have 


(5.146) 


where  it  will  be  assumed  that  the  are  independent.  Substitu- 
ting (5.146)  into  (5.139)  and  averaging  gives 


5-73 


(5.147) 


{ 2*>  2 
Since  EjB  f does  not  necessarily  equal  B , then  in  general  the 

irreducible  error  rate  estimate  is  biased.  The  variance  of  the 

estimate  can  be  easily  found  to  be  given  by 


(5.148) 


where  higher  than  second-order  moments  of  a,  were  considered  to 
be  negligible. 

A comparison  of  (5.145)  with  (5.148)  can  be  made  for  the 
case  where  all  the  B,  are  equal  as  are  all  the  For  this 

case,  (5.148)  reduces  to 


r 

E{Bk 

irr 


} Ek) 


(5.149) 


Since  one  would  expect  that 


(5.150) 


then  for  this  case  the  variance  given  by  (5.145)  will  be  equal  to 
that  given  by  (5.148). 

2 

Using  the  differentiation  technique  to  estimate  B and  for 
nondiversity  operation,  the  fractional  bias  and  estimation  SNR 
can  be  found  from  Section  4.3.2“  to  be  given  by 


\\fhen  the  bias  uncorrected  estimates  of  |T(f,t)|^  are  used  to 
estimate  B^,  then  Appendix  A should  be  used  to  modify  (5.151). 
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(5.151) 


where  Cp  is  the  rms  fractional  error  in  estimating  the  magnitude 
squared  of  the  channel  transfer  function,  K is  the  number  of  inde- 
pendent samples  of  jT(f,t)|2}  and  At  is  the  spacing  between 
samples  of  |T(f,t)p  used  to  approximate  the  derivative.  As 
mentioned  in  Seccion  4.3.2,  we  must  have  tt^B^ (At ) ^ « 1 . 

Figure  5.21  presents  the  performance  of  the  irreducible 
error  rate  estimator  for  first-order  diversity.  In  this  figure, 
it  was  assumed  that  an  independent  sample  of  |T(f,t)|2  can  be 
taken  approximately  every  1/B  seconds;  that  is,  that  K = BT'e, 
where  Te  is  the  time  interval  over  which  the  irreducible  error 
rate  is  estimated. 

5.4.2  Irreducible  Error  Rate  for  FDM-FM  Transmission 

In  this  section  we  will  address  the  problem  of  predicting 
the  irreducible  error  rate  for  FDM-FM  troposcatter  systems.  The 
frequency-selective  fading  causes  the  appearance  of  an  intermodu- 
lation distortion  noise  at  the  discriminator  output  and  will  pro- 
duce an  error  rate  that  cannote  be  reduced  by  increasing  the 
signal-to-noise  ratio,  i.e.,  it  produces  an  irreducible  error  rate. 

Theories  concerning  intermodulation  distortion  in  FM  tropo- 
scatter systems  have  been  developed  independently  by  several 
authors  [5.5]  » [5.7].  We  shall  utilize  the  theory  developed  in 
[5.5],  which  includes  the  effects  of  diversity  combining  and  is 
applicable  to  scatter  communication  channels  with  general  delay 
power  spectrum  Q(£).  As  with  the  others,  this  theory  is  limited 
to  the  case  in  which  the  degree  of  frequency  selectivity  in  the 
rf  signal  bandwidth  is  small. 

The  intermodulation  distortion  is  a realization  of  a sta- 
tionary stochastic  process.  As  such,  it  can  be  described  in 
terms  of  its  probability  density  functions  of  all  orders.  In  par- 
ticular, it  can  be  shown  [5.5]  that  the  function  f(R)  defined  by: 
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(5.152) 


f(R) 


A I 

where  ^3 /AJ  and  n^/&*  are  the  skewness  and  excess  of  the  tropo- 
scatter  channel's  Delay  Power  Spectrum,  exists  such  that,  for 
order  maximal  ratio  predetection  diversity  combining  y( t), 
the  normalized  (unity  rms  value)  second-order  intermodulation 
distortion  has  a first-order  probability  density  function,  W(y), 
given  by: 


w(y) 


MdR 


f 


■/OX)  1 

. 


1 + R2  + 


4f(R) 


-|M+1 

U 


(5.153) 


Equation  (5.153)  for  the  intermodulation  noise  probability 
density  function  when  maximal  ratio  predetection  combining  is  used 
has  been  derived  previously  by  Bello  and  Nelin  [5.5]  and  will  not 
be  repeated  here.  The  probability  density  function  and  cumulative 
distribution  for  selection  diversity  combining  has  been  derived 
in  [5.2],  where  it  was  noted  that  for  M independent  random  vari- 
ables having  a continuous  cumulative  distribution  function  P(y) 
and  a probability  density  function  w(y),  the  probability  that  the 
smallest  of  these  variables,  y^,  is  less  than  a value  T is: 


Prob^  s r}  » 1 - [l  - P(r)]M  (5.154) 

and  the  probability  density  function  associated  with  (5.154)  is: 

wx(y)  = Jnji  -[i-p(r)]Mj^=  m[i  - p(y)]M“1w(y)  (5.155) 

Equations  (5.154)  and  (5.155)  can  be  applied  directly  to 
determine  the  distribution  of  intermodulation  noise  in  a selec- 
tion diversity  FM  system.  The  probability  density  function  cf 
y for  one  branch  is  given  by  (5.153)  with  M=l,  and,  by  inte- 
grating it,  P(y)  can  be  found.  Then,  direct  substitution  into 
(5.154)  and  (5.155)  results  in  the  cumulative  distribution  and 
density  function  of  the  smallest  of  N intermodulation  noises. 

As  a selection  diversity  combiner  is  designed  to  choose  the 
receiver  with  the  smallest  of  the  input  noises,  the  resulting 
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statistics  are  those  of  the  selection  diversity  intermodulation 
noise. 

Equations  (5.153)  and  (5.154)  have  been  evaluated  numeri- 
cally in  {5.2]  and  the  resulting  noise  distributions  are  shown 
in  Figure  5.22.  It  is  seen  that  as  the  order  of  diversity  in- 
creases , the  selection  diversity  cumulative  distribution  lies 
uniformly  to  the  left  of  the  maximal  ratio  cumulative  distribu- 
tion, which  implies  that  the  irreducible  error  due  to  intermodu- 
lation distortion  will  be  lower  for  selection  combining  than  for 
maximal  ratio  combining. 

A quantity  of  importance  in  practical  application  is  the 
ratio  of  the  intermodulation  distortion  power  to  signal  power  in 
a narrow  frequency  band  located  at  some  specified  baseband 
frequency,  f.  Denoting  this  quantity  by  we  may  write 

[5.5] : 

p .(f)  L4y(t) 

- xxp  ,£v  4 (5.156) 

x'  J 

where  L is  tile  rms  multipath  spread  of  the  channel  and,  for  a 
modulating  signal  x(t),  Px^(f)  is  the  power  spectrum  of  x(t)x(t), 
and  P (f)  is  the  power  spectrum  of  x(t).  As  an  example  [5.5], 
consider  the  case  of  a Gaussian  x(t)  with  a flat  power  spectrum 
extending  from  aW  to  W Hz,  and  zero  elsewhere.  For  this  case, 
?7(f,t)  for  aW<f<W  is  given  by 


2(1  - a) 


t2  2-  , 
L rr  f , f 
dev  , 


)2  [MA-W-  4 


(5.157) 


where 


11  , if  aW < f < (l-a)W 
0 , otherwise 


j 1 , if  2aW<  f < W 
1.0,  o therwis  e 


5-78 


Equation  (5.157)  can  be  expressed  as 


- y(t)L4D(f) 


(5.158) 


where 


Thus,  the  rms  multipath  spread  is  the  only  channel  para- 
meter required  to  specify  the  distribution  of  7j(f,t). 

The  error  probability  due  to  intermodulation  noise  can  be 
expressed  in  terms  of  the  conditional  probability  of  error  as 

Pe=/Pefon T)W>dy  <5'159> 

0 

where  W(y)  is  the  probability  density  function  of  the  normalized 
intermodulation  distortion  and  P (•)  is  the  conditional  proba- 
bility of  error.  e 

Figure  5.23  is  a block  diagram  of  a technique  for  estimating 
the  irreducible  error  rate  due  to  frequency  selective  fading. 

The  estimator  uses  estimates  of  the  rms  multipath  spread  to  esti- 
mate the  irreducible  error  rate. 

As  an  example  to  illustrate  the  performance  of  the  estimator, 
we  will  consider  the  systems  analyzed  in  [5.2],  The  systems  are 
called  the  Bell  301B  and  the  Lenlcurt  26B , and  the  irreducible 
error  rate  was  found  in  [5.2]  to  be  given  by  Figure  5.24. 

From  Figure  5.24  the  irreducible  error  rate  can  be  estimated 

by  ^ c 

P.  = A L2  (5.160) 

xrr 

where  A and  C are  constants.  From  Figure  5.24  we  note  that  C 
approximately  equals  the  order  of  diversity.  In  a manner  similar 
to  the  procedure  used  to  evaluate  the  irreducible  error  rate  due 
to  time  selective  fading  [Eqs.  (5.140)  to  (5.145)],  it  follows 
that 
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(5.161) 


2 

where  E|L  j and  o\^  are  the  mean  and  variance  of  the  estimates  of 

■"*2  ^ 

L (see  Section  4.3  and  Appendix  A). 

5.5  F DM/ FM  Error  Ra te  Est imation 


5.5.1  Introduction 

This  section  considers  the  problem  of  estimating  the  error 
rate  of  an  FDM/FM  system.  The  error  rate  for  an  FDM/FM  system 
can  be  approximately  determined  by  separately  finding  the  error 
rate  due  to  nonsolective  fading  (i.e.,  due  to  additive  noise)  and 
the  error  rate  due  to  selective  fading  (i.e.,  due  to  intermodula- 
tion noise).  The  latter  constitutes  an  irreducible  error  rate 
from  the  point  of  view  that  it  cannot  be  reduced  by  increasing 
signal  power.  In  order  to  approximate  the  error  rate  due  to 
additive  noise  and  selective  fading  acting  together,  one  may  sum 
the  error  rates  due  to  these  disturbances  acting  separately. 


5.5.2  Flat  Fading  Error  Rate  Estimation 

In  this  section  we  present  techniques  for  estimating  the 
flat  fading  error  rate  of  an  FDM/FM  system  from  the  received  data 
sigua 1 . 

Using  the  results  of  Rice  [5.8],  Bello  [5.9]  was  able  to 
derive  an  expression  for  the  discriminator  output.  SNR.  He  showed 
that  if  the  typical  signaling  element  is  a sinusoidal  burst  and 
the  instantaneous  predetection  combined  input  SNR,  p,  is  greater 
than  about  1 dB,  the  discriminator  output  instantaneous  SNR  after 
coherent  matched  filter  detection,  s,  is  given  by 


s(p) 


(T 


2(l  - c'p)‘ 


2W(1  - a )e 
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rrp 
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(5.162) 
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where 


a = mean-squared  frequency  deviation  of  modulation 
signal 

B = IF  bandwidth 

r “ B / JVT  (baseband  assumed  flat) 

f , “ location  of  data  subcarrier 
d 

and  the  baseband  extends  from  aW  to  W Hz.  In  calculating  (5.162) 
it  was  assumed  that  there  are  enough  baseband  FDM  signals  so  that 
their  sum  has  close  to  Gaussian  statistics. 

The  SNR  at  the  input  to  the  discriminator  is  distributed 
according  to 


M-l 


W M(p)  « 


M 


(M-l)!fTrk 

k=>l 


for  maximal-ratio  combining  and 

m M“1 

W (p)  ^ 

SVP'  M 

rf  rk 

k=l 


P « rfc  (5.163) 


p « rk  (5.164) 


for  selectioxi  combining.  The  average  error  rate  can  be  found  by 
averaging  the  conditional  error  rate  over  the  distribution  of  the 
instantaneous  SNR.  That  is,  we  can  write 

CO 

P - [ p raCp)]  w(p)  dp  (5.165) 

o 

where  P (s)  is  the  conditional  error  rate, 
e 

Since  the  case  of  p«F  is  the  region  where  most  errors 
occur,  then  the  error  rate  can  be  accurately  estimated  from  the 
distribution  of  the  errors  for  low  SNR.  Therefore,  using  (5.163) 
in  (5.165)  we  obtain 
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p « 
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(5.166) 


j Pefs(p)]  pM_1  dp 
0 


where  A is  a constant  dependent  upon  the  diversity-combining 
technique  and  order  or  diversity.  Therefore,  if  P fs(p)]  is 
known  in  a functional  form,  the  above  integral  can^be  estimated 
and  the  error  rate  can  be  expressed  as 
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(5.167) 


This  expression  has  the  same  form  as  the  error  rates  esti- 
mated in  Section  5.3  [i.e.,  see  Eq.  (5.101)].  Therefore,  the 
error  rate  estimation  technique  proposed  and  analyzed  in  Section 
5.3  can  be  used  to  estimate  the  error  rate  for  FDM/FM  systems 
with  no  intermodulation  distortion.  In  fact,  any  of  the  tech- 
niques analyzed  in  Sections  5. 1.5. 2 and  5.3  can  be  used  to 
estimate  the  error  rate  due  to  noise  alone. 

Several  FDM/FM  systems  were  analyzed  in  [5.2].  The 
AN/GSC-4  is  a 4-phase,  2400  b/s , DPSK  modem.  In  [5.2],  it  was 
found  that  the  conditional  error  probability  for  the  AN-GSC-4  is 


Pefs(p)] 


erfc 


3.33  s(p) 

. 2 \/~6.5  + 3.33  s(p) 


(5.168) 


where  s (p)  is  given  by  (5.162). 

The  Bell  301B  is  a 40.8  kb/s  digital  data  set  that  accepts 
binary  information  at  40.8  kb/s  and  encodes  it  into  4-phase 
DPSK.  The  conditional  error  probability  for  the  Bell  301B  is 
[5.2]  : 


PeU(p)] 
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2.35  s(p) 

_ 2 Vo.5  + 2.35s(p). 


(5.169) 
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Comparing  the  conditional  error  rates  of  the  AN/GSC-4  and  the 
Bell  301B,  one  would  expect  that  the  two  systems  would  have 
virtually  the  same  error  rate. 

The  third  modem  to  be  considered  in  the  Lcnkurt  26B.  The 
Lenkurt  26B  uses  a duobinary  modulation  technique  where  the  in- 
formation is  transmitted  by  an  FSK  modulation  that  deviates 
+ 600  Hz  relative  to  the  center  frequency  of  1.7  kHz.  The  con- 
ditional error  probability  for  the  Lenkurt  26B  is  [5.2]: 


Pc[s(p)] 


3 

4 


erfc 


Vs (p) 


(5.170) 


Thus,  for  the  above  three  FDM/FM  modems,  the  error  rate  due 
to  flat  fading  can  be  estimated  using  any  of  the  techniques  ana- 
lyzed in  Sections  5.1 -5.3. 

5.5.3  FDM/FM  Error  Rate  Estimation  Due  to  Noise  and 
Intermodulation  Distortion 

To  estimate  the  error  rate  due  to  noise  and  intermodulation 
distortion,  one  may  sum  the  flat  fading  error  rate  estimate  of 
Section  5.5.2  with  the  irreducible  error  rate  estimate  given  in 
Section  5.4.2.  That  is,  the  error  rate  for  FDM/FM  transmission 
can  be  expressed  as 


p « p + P 

e irr  FF 


(5.171) 


where  rr  is  the  irreducible  error  rate  due  to  intermodulation 


distortion  and  for  the  Bell  301B  and  the  Lenkui't  26B  systems  can 
be  estimated  by  [see  Eq.  (5.160)] 


P.  - A L 
irr 


(5.172) 


where  A and  C are  constants.  The  flat  fading  error  rate  is 
given  by  (5.167)  and  can  be  estimated  by 


k53! 


(5.173) 
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where  M is  the  order  of  diversity,  is  an  estimate  of  the  mean 
SNR  on  the  lct^  diversity  branch,  andKB  is  a constant  dependent 
upon  the  form  of  the  conditional  error  rate. 

From  (5.171)  - (5.173)  the  error  rate  for  FDM/FM  systems  can 
be  estimated  by 
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(5.174) 


This  estimate  can  be  implemented  from  the  rms  multipath 
spread  and  the  mean  SNR  estimates  proposed  and  analyzed  in 
Section  4.3.  Frequently,  either  the  flat  fading  error  rate  or 
the  irreducible  error  rate  will  dominate  and,  for  these  cases, 
the^convergence  of  Pe  can  be  directly  related  to  the  convergence 
of  Fpp  or  P^rr.  The  estimation  of  Ppp  was  considered  in 
Sections  5.1  to  5.3,  while  Section  5.4.2  considered  the  estima- 
tion of  P. 

irr . 

5 . 6 Estimation  of  Error  Rate  for  High-Speed  Data  Transmission 

5.6.1  Introduction 

In  this  section  we  address  the  problem  of  estimating  error 
rates  for  high-speed  data  transmission  systems  employing  a 
multipath-resisting  receiver  (e.g.,  maximum  likelihood  or  non- 
linear equalizers).  For  these  systems,  in-band  diversity  in- 
creases the  effective  order  of  diversity  and  results  in  an 
improvement  in  performance. 

High-speed  data  transmission  involves  reception  of  signaling 
elements  corrupted  by  intersymbol  interference.  Optimum  demodu- 
lation techniques  exist  (e.g.,  maximum  likelihood  and  nonlinear 
equalizers)  which  are  claimed  [5 . 10]  , [5 . 1.1]  to  come  close  to  eli- 
minating the  effect  of  intersymbol  interference.  The  limiting 
performance  obtainable  is  just  that  associated  with  the  matched 
filter  receiver  for  a single  pulse.  In  this  section  we  discuss 
the  estimation  of  error  probabilities  for  this  ideal  multipath- 
resisting  receiver.  For  practical  receivers,  one  must  assign 
some  degradation  in  SNR.  However,  it  appears  that  this  degrada- 
tion need  not  be  more  than  1 or  2 dB.  Since  the  number  of  prac- 
tical implementations  possible  is  very  large,  it  appears  most 
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meaningful  to  study  the  estimation  of  the  performance  of  the 
optimum  receiver  and  then  apply  correction  factors  for  the  known 
nonideal  performance  of  the  actual  implemented  receiver. 

5.6.2  Single-Pulse  Matched  Filter  System 

The  optimum  performance  we  wish  to  compute  is  associated  with 
a system  in  which  the  receiver  is  a coherent  matched  filter  re- 
ceiver, matched  to  the  response  of  each  pulse  and  completely 
eliminating  the  effect  of  intersymbol  interference.  For  simpli- 
city of  notation,  consider  the  pulse  transmitted  at  t = 0,  and 
let  Ag(t),AG(f)  denote  the  transmitted  pulse  and  its  spectrum, 
respectively.  The  information  is  contained  in  the  value  of  A 
which  takes  on  the  values  +1  according  to  a binary  alphabet. 
Figure  5.25  gives  the  signal  processing  operations  involved.  The 
pulse  duration  is  so  much  smaller  than  the  fading  time  constant 
that  the  channel  may  be  regarded  as  frozen  in  computing  the  pulse 
response.  Confining  our  attention  to  a single  diversity  channel 
for  a moment,  the  spectrum  of  the  received  pulse  is  AG(f)T(f,0) 
and  the  transfer  function  of  the  matched  filter  is  then 
K*(f ) T*(f ,0) . The  complex  sampled  matched  filter  output  is 
given  by 


W - A f |G(f)l2|T(f,0)|2df  + f N(f)  G*(f ) T*(f,0)  df 

(5.175) 


where  N(f'  is  the  complex  spectrum  of  the  white  additive  complex 
Gaussian  noise.  This  spectrum  is  also  white  [5.12], 


N*(f)N(f  + F)  = 2Nq  6 (F)  (5.176) 

where  Nq  is  the  one-sided  power  density  of  the  real  additive 
noise . 

Expression  (5.175)  is  in  the  form 

W = A S + N (5.177) 


where  N is  a complex  Gaussian  noise  whose  strength  is  readily 
computed  to  be 
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Figure  5.25  Single-Pulse  Matched  Filter  System 


(5.178) 


In! 2 - 2Nqs 


where 


S ~ f I G(£) I 2 | T(f , 0)  | 2 df  (..  .179) 

Of  course,  S is  a random  variable  dependent  upon  the  time  of 
transmission,  but  the  fading  is  so  slow  relative  to  the  bit  rate 
that,  for  purposes  of  analysis,  we  may  consider  S fixed  for  many 
bits.  Thus,  we  may  evaluate  bit  error  probabilities  by  assuming 
S fixed  and  then  subsequently  averaging  over  the  statistics  of  S. 

When  diversity  channels  with  independent  additive  noises 
are  included  in  the  analysis,  it  is  readily  seen  that  we  must 
redefine  S as 


9 M 9 

s-  /lC(E)r  £ iT.  <£,0)1 ‘‘df  (5.180) 

1 k=l  K 

th 

where  Tj<.(f,0)  is  the  transfer  function  of  the  k diversity  chan- 
nel and  M is  the  order  of  diversity.  Equation  (5.178)  still 
applies  but  with  the  new  definition  of  S in  (5.180). 

The  Re{  } operation  in  Figure  5.25  corresponds  to  the  co- 
herent detection  operation  of  the  PSK  system.  Thus,  prior  to  the 
binary  (^  or  >)  decision,  we  use  the  variable 


U = Re{W]  = AS  + Re{N] 


(5.181) 


The  variable  Re{N]  is  real  Gaussian  with  variance 


Re2{N]  = ~ I N I 2 


(5.182) 


5.6.3  Estimation  of  S(t) 

In  order  to  estimate  the  bit  error  probability,  the  estima- 
tion of  S(t)  must  be  considered,  where  S(t)  is  given  by 

M 

s(t)  = f 1 g (f ) | 2 It  (f,t)|2df 
J k=l  K 
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(5.183) 


Two  techniques  for  estimating  S(t)  from  the  received  data 
signal  may  be  used.  The  first  approach  uses  samples  in  frequency 
of  an  estimate  of  the  magnitude  squared  of  the  channel  transfer 
function  to  express  the  integral  as  a finite  sum.  This  approach 
is  presented  in  Figure  5.26.  Therefore,  the  estimate  of  S(t  ) is 
given  by  p 
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G(Fq) I 


(5.184) 


where  W = Q • £F  is  the  bandwidth  of  G(f). 

With  the  estimate  of  the  magnitude  squared  of  the  channel 
transfer  function  given  by 


!Tk(Fq»tp)!  = !T(F„,tJ|Z  (!  + €,_)  + 6 


q p 


kq' 


kq 


(5.185) 

where  6^,  ^q  are  zero  mean  and  independent  estimation  errors 
with  (see  Section  4.1) 


2 2 

m c (5.186) 

€kq 


where  e and  cg  are  given  by  (4.57)  and  (4.100),  respectively. 
From  (5.184)  and  (5.185),  the  mean  and  variance  of  S(tp)  condi- 
tioned upon  the  channel  transfer  function  are  given  by 
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Figure  5.26  Estimation  of  Total  Received  Signal  Power 


For  large  Q,  the  summation  of  (5.187)  converges  to  the  de- 
sired integral  given  by  (5.183)  and  the  variance  of  S(t  ) becomes 
small.  Thus,  when  Q is  large,  a good  estimate  of  S(tp)^can  be 
obtained.  However,  a large  Q means  that  many  magnitude- squared 
channel  transfer  function  estimators  would  be  required. 


The  second  approach  for  estimating  S(t)  that  will  be  con- 
sidered utilizes  the  technique  employed  to  estimate  the  magnitude 
squared  of  the  channel  transfer  function.  The  estimator  to  be 
considered  is  given  by  Figure  5.27  and  has  the  same  functional 
representation  as  the  magnitude- squared  channel  transfer  function 
estimator  analyzed  in  Section  4.1.  The  averaging  filter  duration 
will  be  chosen  long  enough  to  average  out  data  fluctuations,  yet 
short  enough  not  to  filter  the  channel  fluctuations.  The  band- 
width of  h(t)  should  be  large  enough  to  pass  the  signal  without 
distortion,  but  not  so  large  that  the  amount  of  noise  power 
passed  is  excessive. 


The  estimate  of  S(t)  averaged  over  the  data  and  noise  can  be 
given  by 


E{s(t)}  = K(0) 


1 H (0 ) { 2 J | G(f  ) 1 2 | T(f  ,t)  | 2 df 


+ 2N0  / |H(f ) 1 2 df 


(5.189) 


where  H(f)  was  assumed  to  be  flat  over  the 


duration  of  G(f)  and 


K(0)  = f k(t ) dt 


(5.190) 


The  bias  and  variance  of  the  estimator  in  Figure  5.27  can  be  com- 
puted as  was  done  for  the  estimation  of  lT(f,t)|2.  Frequency- 
selective  distortion  can  be  neglected  because  an  integral  of 
|G(f)|2|T(f,t) j2  over  f is  computed.  Thus,  the  output  averaging 
filter  can  be  adjusted  to  minimize  the  effects  of  time-selective 
fading.  It  appears  that  the  desired  time-variant  energy  can  be 
estimated  more  accurately  than  |T(f,t)|2  so  that  the  second 
approach  is  preferred. 
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Figure  5.27  Estimation  of  Received  Signal  Power 


Once  S(t)  is  estimated  by  either  of  the  two  methods  presented 
above,  the  error  rate  may  be  estimated  from  the  conditional  bit 
error  probability.  Figure  5.28  presents  a technique  for  esti- 
mating error  rate  from  estimates  of  S(t).  As  was  shown  in 
Section  5.1.4,  the  conditional  error  rate  can  also  be  used  to 
estimate  short-term  error  bursts. 

5.7  Error  Rate  Estimation  Based  Upon  Interference 

In  this  section  we  describe  the  way  that  error  rate  estima- 
tion is  affected  by  errors  in  estimating  the  interference  para- 
meters and  channel  transfer  function.  The  error  rate  predictions 
used  here  are  in  the  nature  of  upper  and  lower  bounds  on  proba- 
bility of  error.  This  general  approach  is  used  here  because  of 
the  wide  variety  of  interference  that  might  occur  on  all  of  the 
channels  discussed  in  this  report.  Interference  could  take  the 
form  of  man-made  noise,  naturally-occurring  interference,  and 
deliberate  jamming.  Due  to  the  unpredictability  of  the  inter- 
ference, it  is  not  feasible  to  attempt  precise  error-rate  predic- 
tions at  a given  SNR,  and  any  results  relating  to  specific 
interference  models,  such  as  normally  distributed  interference, 
would  be  of  extremely  limited  utility. 

In  the  work  reported  below,  the  additive  interference  will 
be  allowed  to  have  general  statistics.  Restrictions  on  the  inter- 
ference will  be  minor.  At  first,  we  will  describe  how  error  rate 
bounds  can  be  developed  by  using  as  the  only  constraints  a speci- 
fication of  the  mean  value,  a,  and  an  upper  bound  on  the  values 
that  the  generalized  noise  variable  can  take  on.  Later,  in  this 
section,  we  will  replace  the  upper  bound  constraint  with  a 
constraint  on  p,  the  probability  that  the  interference  variable  U 
exceeds  some  fixed  threshold  U* . The  motivation  for  this  new 
approach  arises  from  the  simplicity  of  the  schemes  for  estimating 
p and  the  reliability  with  which  it  can  be  measured,  as  discussed 
earlier  in  Section  4. 5. 4. 3.  In  the  development  leading  to  both 
of  these  schemes,  the  fading,  as  before,  is  assumed  to  be  charac- 
terizable  as  a comp lex- valued  Gaussian  process  and  is  assumed  to 
be  sufficiently  slow  and  to  have  a sufficiently  small  degree  of 
frequency  selectivity  so  that  a received  signaling  element  is 
undistorted  by  the  channel  and  differs  from  the  transmitted 
signaling  element  only  in  a random  amplitude  change  and  phase 
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Figure  5.28  Error  Rate  Estimator  for  Multipath  Resisting  Receiver 


shift.  System  degradations  due  to  fading  that  is  not  restricted 
in  this  way  have  been  treated  elsewhere  in  this  report.  Since 
the  effect  of  these  distortions  can  be  reduced  to  the  point  where 
system  degradation  is  determined  primarily  by  nonseiective  fading 
and  additive  interference,  it  is  this  last  situation  that  we  con- 
centrate on.  Hence,  with  s(t)  denoting  the  complex  envelope 
representation  of  a signaling  element  of  duration  T received  in 
the  time  interval  0<t<T  for  a unit  gain  channel  and  no  additive 
noise,  we  have  for  the  received  signaling  element  r(t),  including 
fading  and  additive  interference,  the  representation 


r(t)  * G s(t)+n(t)  0<t<T 


(5.191) 


where  G is  a complex- valued  normally  distributed  variable,  and 
n(t)  is  the  additive  interference. 

It  was  pointed  out  in  Sec tion  4 . 5 that  three  important  matched 
filter  receivers  are  characterized  well  by  detector  operation 
that  can  be  represented  as  the  comparison  of  a quadratic  form 
with  a zero  threshold.  The  quantity  q^  is  given  by 
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(5.192) 


where  L is  the  order  of  diversity,  G is  the  complex- valued 
normally  distributed  signal  term  on  the  pfch  diversity  branch, 
and  £p  and  77 p are  complex  noise  terms  due  to  the  p^  diversity 
receiver.  Specific  forms  were  given  in  Table  4.5.  For  diversity 
channels  that  fade  independently,  but  with  identical  statistics, 
the  random  variables  jGpj-  become  mutually  independent  and  identi- 
cally distributed.  r 

For  future  use,  we  note  no  loss  of  generality  accrues  if  the 
normalization 


(5.193) 


is  used  which  means,  effectively,  that  the  quantities  U and  V 
are  "instantaneous"  noise  to  average  signal  power  ratios. 
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As  shown  in  [5.17],  the  problem  of  determining  error  rate 
bounds  is  a double  integral  extremum  problem  that  involves  the 
interference  only  through  the  joint  probability  density  function 
of  the  two  interference  variables,  U and  V, 


- w 

P=1 

(5.194) 

L 

’ - S ' v* 

(5.195) 

The  upper  bound  on  the  double  integral  is  achieved  when  U and  V 
are  assumed  to  have  nonzero  values  mutually  exclusively,  i.e., 


UV  - 0 


(5.196) 


and  the  lower  bound  is  achieved  when 


U 


(5.197) 


It  is  demonstrated  in  [5.17]  that  the  error  probability  is 
bounded  in  the  following  way: 


/F  (A)  Gt  (A)  dAsp  s /F  (A)  H (A)  dA 


(5.198) 


where  Fj(A)  is  the  probability  density  function  for  A,  the  arith- 
metic mean  of  U and  V. 
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are  given,  respectively,  by 
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and 


GL(A) 


L-l 

e"A  £ I (A) 
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(5.201) 


where  I^(A)  is  the  modified  Bessel  function  of  order  r. 

For  the  general  situation  being  treated  here,  the  pdf  Ft  (A) 
is  not  known.  It  is  clear,  though,  that  a maximization  of  the 
right-hand  side  of  (5 .198)  and  a minimization  of  the  left-hand  side 
leads  to  reliable  bounds  on  error  probability,*  This  extremization 
problem  only  requires  proper  choice  of  the  minimizing  and  maximi- 
zing pdf's,  a task  which  becomes  extremely  simple  through  applica- 
tion of  the  "Completeness  Theorem"  proved  in  [5.17],  A reasonable 
approach  to  the  problem,  of  course,  requires  the  application  of 
constraints,  without  which  the  bounds  would  be  useless.  The 
physical  basis  for  all  of  the  discussion  in  this  section  is  that 
the  constraints  can  be  generated  by  direct  interference 
measurement . 


5.7.1  Bounds  on  Error  Rate 

In  order  to  plot  the  bounds  on  error  rate  as  a function  of 
signal-to-noise  ratio,  it  is  necessary,  here,  to  first  define  a 
composite  signal-to-noise  ratio. 
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where  the  normalization  in  (5.193)  has  been  used  and  is  defined 
as  the  mean  of  the  noise  variable  U or  A,  i.e.. 
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(5.203) 


We  point  out  that  coherent  PSK  systems,  for  which  U=V,  require 
only  that  the  left-hand  side  of  (5.198)  be  dealt  with,  i.e., 
minimize  and  maximize  J*  F^(A)  G^(A)  dA  to  get  both  upper  and 
lower  bounds  on  the  error  rate. 
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where  A is  just  the  arithmetic  mean  of  U and  V: 


U+V 
A - — 


There  are  two  constraints  that  are  used  in  establishing  the 
error  bound.  The  first  of  these  is  the  maximum  allowable  peak 
noise  power  to  average  signal  power  constraint  defined  below s 


y = Max 


Peak  U 
I 2L 


; U = V 


; UV  - 0 


(5.204) 


The  second  is  the  maximum  allowable  peak-to-average  interference 
constraint  or  "crest  factor": 
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(5.205) 


In  both  of  these  expressions.  Peak  A refers  to  the  peak  or  maxi- 
mum value  of  the  variable  A.  The  first  thing  to  be  noticed  about 
these  expressions  is  that  y^  does  not  depend  on  a [or  signal-to- 
noise  ratio;  see  (5.202)]  whereas  /3^  doe  . This  is  because  the 
earlier  constraint  arises  because  of  . ae  finite  dynamic  ranges  of 
receivers,  which  cause  the  outputs  to  saturate  for  sufficiently 
high  input  voltages,  and  it  is  safe  to  assume  that  detected 
noises  cannot  exceed  the  design  value  of  output  signal  power  by 
more  than  some  factor,  i.e.,  there  is  'me  factor  y^  for  which 
Peak  U is  always  less  than  2>'lL.  This  factor,  alone,  places 
limits  on  the  integrals  in  (5 . 198 ) which  cannot  be  exceeded;  this, 
in  itself,  is  enough  to  give  meaning  to  the  extremization  problem. 
Unfortunately,  the  error  bounds  that  would  result  from  such  a 
formulation  would  be  extremely  weak. 

The  use  of  fir  has  the  effect  of  further  constraining  Peak  U, 
and  has  the  end  effect,  typically,  of  tightening  the  bounds  on 
probability  of  error.  As  opposed  to  the  constraint,  yL,  which 
arises  mainly  from  consideration  of  receiver  limitations,  the 
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constraint  /?l  arises  solely  from  consideration  of  the  interfer- 
ence statistics.  The  use  of  some  specific  value  requires 
some  a priori  knowledge  of  these  statistics.  This  knowledge 
could  be  obtained  by  utilizing  realistic  models  of  naturally 
occurring  or  man-made  interference.  Unfortunately,  the  interfer- 
ence is  often  of  such  an  unpredictable  nature  that  the  only 
reliable  way  to  obtain  any  useful  information  about  interference 
statistics  is  through  measurement.  This  fact,  of  course,  was 
the  motivation  for  the  measurement  schemes  discussed  in  Section 
4.5  and  much  of  the  work  to  follow. 

We  now  review  the  results  of  the  extremization  procedures 
as  reported  in  [5.17]. 

First  we  recall,  as  discussed  earlier,  that  the  left-hand 
side  of  (5 . 198 ) must  be  minimized  and  the  right-hand  side  maximized. 
In  the  case  U=V,  we  need  only  deal  with  the  left-hand  side  of 
(5.198).  This  gives  us  the  three  extremization  problems  listed 
below: 


I.  (UV=0)  Maximize  the  right-hand  side  of  (5.198)  to 
get  pJ^ax. 

II.  (U=V)  Maximize  the  left-hand  side  of  (5.198)  to 
get  p™ax. 

Jj 

III.  (U=V)  Minimize  the  left-hand  side  of  (5.198)  to 

get  p™in. 

L 

The  maximization  problem  for  the  UV=0  case  is  to  maximize 

A 

al 

PL  “/  FL(A)HL(A)dA  (5.206) 

0 


subject  to  the 


a. 


L 


constraint 


/ AFL(A)dA 
0 


where 


A 


Min 


(5.207) 


(5.208) 


5-101 


For  the  U-V  case  (problems  II  and  III  above)  the  extremization 
problem  involves 


pL  “ / WL(U)GL(U)dU  (5.209) 

0 


subject  to  the  constraint 


«L  » j U WL(U)  dU 

0 


where 


UL  Min|  2L  rL  , 2Al 


(5.210) 


(5.211) 


The  results  of  the  extremizations  as  a function  of  p,  the 
SNR,  are  given  by 


Max 


(2//Sl)  • Hl(L  0L/p)  ; 
(2L/p)  • HL(pL)/JuL  ; 
HL(2L/p) 


PiL  R-J  Ml 

( 2L/tu  j)  < p £ L B^J  i il 
2 /y^sps  2L/aL 


(UV=0) 


(5.212) 


where  is  the  point  of  tangency  between  a straight  line  emana- 
ting from  the  origin  and  the  Hl(A)  curve,  and  has  the  values 
M-^O,  ^2=51*®»  an<3  ji^-4.85.  For  problem  II,  the  results  are: 


max 


(1/«L)  • GL(2L|3L/p);  ps  2L/3l/i;l 
(2L/p)  • ’ 2L/z^L  <:  p ^ 2L  j3T /i^L 


GL(2L/p) 


l/rL^  p <;  2L/ 


(U-V) 

(5.213) 
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who 


re  is  analogous  to  /j^;  i.e., 


C]SV\),VL  " Max{GL(U)/u} 


(5.214) 


and  has  the  values  ^“l*^,  p^=S.S. 

The  minimization  in  problem  III  is  slightly  more  complicated 
and  requires  definition  of  a^,  where  cr^  is  such  that 


w. : w 

0L  ’ °L 


= Max 
U, 


w : w 

\ “ UL  J 


(5.215) 


Since  crT  is  a function  of  U , it  is  convenient  to  define 
L 


"l  “ W 


(5.216) 


The  lower  bound  on  error  probability,  as  a function  of  SNR,  is 
given  by 


min 


GL<2Ln>-  2L(’'L-?)SL(2L’'L); 


p < VrL 


dL(2L/S L/p)  - (2L/o)03l-  l)SL(2LjSL/p); 


GL(2L/p) 


j8L/yL <P<  2L/01 
p > 2L/0t 


L > 1 

L > 1 
L > 1 


(5.217) 


where  0T  is  the  solution  to  the  equation 
L 


x = fL(PLx) 


(5.218) 
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ttv- ■ rr  mi  .if..  w-ui 


and  where  we  have  defined  the  maximum  slope  as 


W 


gA>  - GJvv 

°L  - flA> 


(5.219) 


All  of  the  results  catalogued  above  take  ca  considerably 
simpler  form  for  the  nondiversity  case,  L=l.  We  have 


Max  N 

Pi  - H.  (2/p) 


P s 2/y1 


max 


G1(2/p) 


p * l hx 


min 


?!  “ (2/p)-G1(2y1)/2y1  /y^ 


= (l/01)*G1(2fl1/p) 


p * /31/y1 


(5.220) 


Seme  numerical  results  for  different  types  of  modulation  and 
for  both  nondiversity  and  dual  diversity  operation  are  shown  in 
Figures  5.29-5.31.  The  details  of  each  physical  situation  are 
entered  on  each  figure. 

5.7.2  Reliability  of  Bounds 

5. 7. 2.1  Preliminary  Comments 

The  curves  presented  in  the  previous  section  indicate  that 
with  an  extremely  limited  knowledge  of  noise  statistics,  i.e., 
the  peak  noise-to-average-signal  ratio,  y^,  and  the  crest  factor, 
one  can  develop  meaningful  bounds  on  probability  of  error  as 
a ^function  of  the  signal-to-noise  ratio.  In  this  section,  we 
are  not  so  much  interested  in  how  well  one  can  generate  curves  of 
this  type  as  we  are  in  the  question  of  determining  how  well, 
given  a set  of  interference  measurements  (and  channel  transfer 
measurements),  we  can  find  just  the  two  numbers  which  bound  the 
error  probability  under  the  true  noise  condition.  A schematic 
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Figure  5.30  Error  Probability  Bounds  for  Bina 
Function  of  Hatched  Filter  Output 
Diversity  with  Square  Law  Combini; 
(From  Ref.  [5.17]) 


of  the  processing  operation  is  given  in  Figure  5.32.  The  inter- 
ference measurements  are  just  those  discussed  earlier  in 
Section  4.5.  The  measurements  of  the  mean  value  of  the  magnitude- 
squared  transfer  functions  can  be  obtained  through  use  of  the 
estimation  schemes  utilizing  received  data  only,  as  discussed 
earlier  in  Section  4.1  or,  alternatively,  can  be  obtained  by 
making  use  of  the  special  probing  signals  discussed  in  Section  6. 
The  exact  details  of  the  manner  in  which  this  latter  measurement 
is  obtained  will  not  be  of  large  importance  in  the  following 
discussion. 

As  illustrated  in  Figure  5.32,  the  estimates  p™ax  and  p™^n 

A A ■■■"  1 m' 

are  determined  by  processing  or,  dr,  and  |G|^.  The  estimates  are 

A -Ll  A. 

used  in  different  ways;  and  | G|?  are  used  to  estimate  the 
signal-to-noise  ratio,  p.  This,  of  course,  specifies  tie  point 
at  which  the  function  forms  in  (5.212)-  (5 . 220)  must  be  calculated. 
The  estimate  of  /3L  enters  as  a parameter  affecting  both  the  value 
of  the  bound  and  also  the  interval  (and  corresponding  functional 
form)  that  must  be  used;  e.g.,  see  the  first  line  of  Eq.  (5.212). 
The  numerical  processing,  of  course,  just  refers  to  numerical 
implementation  of  the  results  presented  in  Section  5.7.1. 

Finally,  we  point  out  that  the  emphasis  in  the  remainder  of 
this  section  is  on  examining  the  feasibility  of  developing  error 
bounds  from  interference  measurements.  It  is  not  necessarily 
true  that  the  analytical  bounds  presented  in  Section  5 . 7 . 1 consti- 
tute the  best  approach  to  the  problem  from  the  viewpoint  of 
utilizing  interference  measurements.  One  reason  for  this  relates 
to  difficulties  associated  with  estimating  the  interference 
maximum  (crest  factor)  as  discussed  in  Section  4.5.  It  would 
seem  that  alternative  approaches  to  the  general  problem  of  bound- 
ing error  rate  could  yield  bounding  curves  that  make  more  effec- 
tive use  of  the  interference  measurements.  An  alternative 
approach  that  does  not  make  use  of  the  crest  factor  constraint 
is  discussed  in  Section  5.7.3. 

5. 7. 2. 2 Bound  Errors  Arising  from  Errors  in  Interference 
Measurement 

We  now  examine  the  way  that  the  error  bounds,  discussed  in 
Section  5.7.1  and  illustrated  in  Figures  5 .29  - 5 .31 , are  affected  by 
errors  in  the  interference  measurements.  As  a first  step  in  this 
direction,  we  concentrate  on  determining  the  error  in  the  signal- 
to-noise  ratio  measurement. 
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The  problem  of  determining  the  degree  of  accuracy  attainable 
in  estimating  the  signal-to-noise  ratio,  p,  is  central  to  evalu- 
ating any  schemes  that  use  interference  measurements  to  develop 
error  rate  bounds.  The  reliability  of  the  signal- to-noise  ratio 
estimate  depends,  of  course,  on  the  accuracy  attained  in  esti- 
mating the  mean  values  of  the  channel  and  interference  parameters 
in  (5.202).  The  assumption  of  identical  statistics  for  G on  each 
diversity  branch  implies  that 


(5.221) 


It  was  the  assumption  of  identical  interference  statistics  on 
each  diversity  branch  that  led  to 


P 


t|lG|2 


O', 


(5.222) 


and  finally  the  normalization 


||G|2  » 1 (5.223) 

that  led  to  the  right-hand  side  of  (5.202). 

Here,  we  note  from  the  discussion  of  generalized  noise  vari- 
ables in  Section  4.5  that 


G = g(0) T (5.224) 

where  the  received  signal  is  given  by 

y(t)  - g(t)  s(t)  + n(t)  (5.225) 

where  g(t)  just  represents  the  time-varying  complex  gain,  i.e., 
the  channel  transfer  function.  Thus, 


T2  |g(t)|2 


(5.226) 


and  its  estimate  is  given  by 


A 

I G|  2 = T2  )g(t>]  2 (1  + cT) 
- I G|  2 (l  + cT) 


(5.22 

(5.227) 


where  just  represents  the  error  in  estimating  the  mean-squared 
value  of  the  channel  transfer  function.  This  error  has  been  dis- 
cussed at  great  length  in  Section  4.1  and  elsewhere  in  this 
report.  Finally,  using  the  normalization  in  (5.223)  and  the  ex- 
pression for  p in  (5.222),  we  obtain  as  our  estimate  for  p 


P 


“l 


I1'  <1  + fT) 

aT 


(5.228) 


It  is  convenient  to  define  fractional  errors  for  p and  ol 
and  to  denote  them  and  ea,  respectively.  These  random  vari- 
ables are  defined  according  to 


P 


p(1  + fl°)  = ^ (1  + fP) 


(5.229) 


and 


a * a(l  + <ra)  (5.230) 


Now,  substituting  (5.230)  in  (5.228)  and  comparing  the  result  with 
(5.229) , we  obtain  the  fractional  error 
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(5.231) 


CT  “ ?a 


P 1 + ^ 


for  the  signal-to-noise  ratio. 

The  other  major  factor  on  which  the  expressions  for  proba- 
bility of  error  bounds  in  (5.212)  and  (5 . 213)  depends  is  the  crest 
factor  0^.  Estimates  of  this  quantity  are  formed  from  estimates 
of  and  Umax,  the  maximum  value  of  U,  as  discussed  in  Section  6. 
In  fact,  we  recall  from  (4.532)  that  B ^ has  the  form 


„ A „ (Peak  U 
B - Max<- — 


(5.232) 


Defining  a normalized  error  Co  for  the  8^  estimate,  we  ha>r 


«L  " \(1  + V 


(5.233) 


and  for  the  normalized  error  in  estimating  U 


IJ  - U (1  + c ) 
max  max  u 


(5.234) 


The  8^  estimate  is  obtained  numerically  by  dividing  the  U esti- 
mate by  the  a estimate,  i.e., 

L 


= 


(5.235) 


Hence,  using  (5.230),  (5.234),  and  (5.235),  we  obtain 


e = 1 a 

* 1 + 


(5.236) 
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Thus,  in  (5.231)  and  (5.236)  we  have,  respectively,  expressions  for 
the  fractional  measurement  errors  in  measuring  the  signal-to-noise 
ratio,  p,  and  the  crest  factor, 

In  order  to  focus  on  a goal  for  the  calculations,  we  will 
concentrate  on  the  reliability  with  which  the  upper  bounds  in 
(5.212)  can  be  determined  from  the  measurements.  We  note  that  the 
results  that  will  be  obtained  cannot  differ  much  from  those  that 
would  follow  from  a similar  analysis  applied  to  (5.213);  in  fact, 
the  fractional  measurement  error  in  determining  p^ax  and  pj*ax  is 
exactly  the  same  for  the  "middle  range"  of  p values.  For  the 
other  ranges,  the  results  would  be  similar  because  of  the  very 
similar  smooth  behavior  of  Hl(x)  and  Gl(x).  For  the  lower  bound, 
the  situation  is  not  as  clearcut.  This  is  due  to  the  complexity 
of  the  expressions  in  (5.218).  We  point  out,  however,  from 
Figures  5.42  - 5.45,  that  the  slope  of  the  lower  bound  is  nowhere 
steeper  than  it  is  for  the  upper  bound.  This  is  a rough  indica- 
tion that  sensitivity  to  errors  in  the  p estimate  will  be  less 
than,  or  comparable  with,  that  for  the  upper  bound.  We  also 
point  out  that  the  absolute  errors,  all  other  things  being  equal, 
will  be  less  for  the  lower  bound  than  the  higher. 

We  now  turn  to  direct  calculation  of  the  error  in  setting 
the  bounds.  We  first  deal  with  the  interval  p^L/Si^/p^.  For 
these  large  p values,  we  make  use  of  the  Taylor's  expansion  of 
H^(x)  about  x = 0: 


Hl(x) 


1 L 1 
2 X L! 


1 L+l  L 

2 X (L+l)! 


(5.237) 


Utilizing  this  approximation,  we  obtain 


P 


Max 

L 


J_ 

L 

P 


^L1 

(L-  1)! 


(5.238) 


Now,  using  the  estimates  in  (5 . 229 ) and  (5.233)  to  form  an  estimate 
of  p^ax  a(-  these  larger  p values,  \e  obtain 
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*Max 
Pt  a ~7 


T L-l^L-1 
L ^ 


L (L  - 1) ! 


f(W 


(5.239) 


where 


fW 


(1  + fa) 

<1  + V' 


L-l 


(5.240) 


A Taylor’s  expansion  of  this  function  about  f = 0 and  e =0 


gives  to  first  order  the  result 


(3 


f(<rfl,€  ) * 1 + <L“De  " L« 

P P 0 P 


(5.241) 


Thus,  the  upper  bound  on  probability  of  error  has  the  form 


AMax  max./,  , N 
Pr  « PT  (l  + *i) 


(5.242) 


where  the  fractional  error  is  given  by 


= <L-1)f«  - LfP 


(5.243) 


A 

With  the  totally  reasonable  assumption  that  the  $ and  p measure- 
ments are  weakly  correlated,  we  obtain  finally  as  the  mean-squared 
squared  fractional  error  the  result 


2 < 
€1  w 


(5.244) 


For  the  "middle  range"  of  p values,  (2L/p^)  ^ p s:  L we 
note  that  the  upper  bound  probability  of  errorTias  the  gener'al 
form 


AMax 
P 


const (2L/p) 


(5.245) 


This  is  also  true  for  the  middle  range  in  (5.213). 
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Writing 


AMax  Max,-  , v 
p - p <1  + c2> 


it  is  a simple  matter  to  determine  that 


(5.246) 


(5.247) 


We  now  turn  to  the  range  of  p consisting  of  the  smallest  values  of 
of  p,  2/7^  ps:  2L/Ml,  we  form  [see  (5.212)]  the  estimate 


'Max 


* ht  ■ 
L / 


2L 


p(1  + V. 


(5.248) 


With  a Taylor's  expansion  with  respect  to  we  obtain 


'Max 


m?) 


dH 


L /2L\ 

dx  \ p / Cp 


(5.249) 


The  largest  value  of  p in  the  interval  of  concern  is  2L/p^.  At 
this  value  of  p,  the  factor 


(5.250) 


takes  on  unity  value  (by  definition  of  Ml);  it  is  smaller  for  all 
other  values  of  p.  It  follows  immediately  from  the  representation 


that 


*Max 
P 


Max,-  , v 

pl  d+«3) 


2 2 

e o « 

3 p 


(5.251) 
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is  a "worst-case"  approximation  for  the  mean-squared  fractional 
error. 

From  (5.231)  and  (5.236)  we  now  utilize  the  reasonable 
approximations 


2 

<6 


& 


(5.252) 

(5.253) 


This  allows  us  to  represent  our  present  results  in  terms  of  the 
results  derived  earlier  in  Section  4.5.  We  obtain 


2 


p 2 L/3L/pL 


(2L/  fi^)  spi  L $ f juL 


(H 


2/yh<  p*  2L/pL  (5.254) 


It  is  clear  that  accurate  measurements  of  the  interference 
mean  and  maximum  values  in  combination  with  accurate  measurements 
of  the  channel  transfer  function  magnitude-squared  value  can  lead 
to  accurate  determinations  of  the  bounding  curves  for  error  rate. 
To  give  a practical  feoL  for  the  effect  of  the  measurement  errors 
on  the  bounding  curve,  we  have  included  in  the  upper  left-hand 
corner  of  Figure  5.33  a plot  of  the  95%  confidence  limit  versus 
rms  fractional  error  in  estimating  the  bounding  probability. 
Also  shown  in  this  figure  are  the  95%  confidence  limits 
for  the  pulse  noise  jamming  signal  discussed  in  Section  4.5.  It 
is  reasonable  to  expect  for  the  naturally  occurring  types  of 
interference,  that  considerable  improvement  could  be  obtained. 

Finally,  we  point  out  that  the  confidence  limits  at  the 
break  points  cannot  be  determined  exactly  by  analytical  means. 

It  is  safe  to  assume,  however,  that  the  largest  of  the  confidence 
limits  on  either  side  of  a juncture  apply  over  a small  region 
about  the  juncture;  the  width  of  this  region  depends,  of  course, 
on  the  accuracy  with  which  the  signal- to-noise  ratio  can  be 
estimated. 
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Figure  5.33  Confidence  Limits  (95%)  for  Error  Probability 
Bounding  Curves 
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5.7.3  An  Alternative  Approach 

5.7. 3.1  Introduction 

Though  estimates  of  8^  and  can  be  used  with  reasonable 
accuracy  to  develop  error  rate  bounds  by  using  the  functional 
relations  in  Section  5 . 7 . 1 , it  is  not  necessarily  true  that  the 
relations  in  Section  5 . 7 . 1 make  the  most  effective  use  of  the  inter- 
ference measurements.  One  of  the  major  difficulties  with  the 
earlier  computational  scheme  was  that  it  required  estimates  of 
the  probability  density  function  maximum.  Depending  on  the 
amount  by  which  the  pdf  tapers  off  at  higher  values,  the  occur- 
rence of  a value  close  to  the  maximum  can  be  a very  low  proba- 
bility event.  This  problem  was  discussed  earlier  in  connection 
with  the  work  of  Section  4.5.  The  use  of  the  measurement  maximum, 
as  a constraint  for  the  extremization  problem,  was  discussed  in 
Section  5.7.1.  Here  we  discuss  the  use  of  a different  constraint. 
Within  a reasonable  range  of  operational  choices,  this  new  con- 
straint relates  to  events  of  higher  probability  than  was  the 
previous  case. 

We  also  note,  for  schemes  aimed  at  measuring  Umax>  that 
there  is  really  no  way  of  knowing  how  accurate  the  measurements 
are  without  some  a priori  knowledge  of  the  probability  density 
function.  Without  some  other  supporting  measurements,  the 
operator  of  the  system  is  to  some  extent  working  in  the  dark;  he 
has  no  knowledge  of  the  quality  of  the  measurements,  nor  how  much 
more  time  would  be  required  to  obtain  more  reliable  measurements. 
This  is  in  contrast  with  the  simple  scheme  for  obtaining  proba- 
bility estimates  discussed  in  Section  4. 5. 4. 3.  We  recall  that 
the  scheme  presented  there  had  as  its  goal  the  measurement  of  the 
probability,  p,  that  the  generalized  noise  variable,  U,  exceeds 
some  fixed  threshold,  U1  . The  rms  fractional  error  was  given  by 


(5-255) 

A 

We  see  that  the  measurement  p itself  provides  some  information 
about  the  number  of  samples  required  to  obtain  a prescribed  ac- 
curacy in  measurement.  This  leaves  the  operator  some  options, 
either  increasing  M or  decreasing  the  threshold,  U1 , thus  increas- 
ing p and  decreasing  the  error  in  (5.255). 

The  extra  flexibility  discussed  above,  in  conjunction  with 
the  superior  convergence  of  the  p estimate  (as  compared  to  the 
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estimate)  motivates  the  use  of  the  p measurement  to  develop 
bounds  on  error  rate  for  various  communication  systems. 

5. 7. 3. 2 Bounds  on  Error  Rate:  Problem  Statement 

As  pointed  out  in  Section  5 . 7 . 1,  the  problem  of  bounding 
probability  of  error,  in  its  most  general  context,  is  framed  in 
terms  of  the  integral 


A 

A, 


-/ 


Fl(A)  Hl(A)  dA 


0 


(5.256) 


which  is  a precise  expression  for  probability  of  error  under  the 
mutually  exclusive  interference  condition,  UV  = 0,  and  in  terms  of 
the  integral 


P 


L 


WL(U)  Gl(U)  du 


(5.257) 


which  is  the  precise  expression  for  probability  of  error  under 
the  identical^ interference  condition,  U = V.  In  both  of  these 
expressions,  Al  and  Ul  are  appropriately  determined  upper  limits 
to  be  discussed  below.  In  the  first  integral,  A is  just  the 
arithmetic  mean  of  U and  V.  In  both  of  the  integrals,  Ht(x)  and 
Gj^x)  are  the  known  functions  given  earlier  in  Eqs . (5.200)  and 
(5 . 201) .respectively , and  Fl(x)  and  Wp(x)  are  the  unknown  proba- 
bility density  functions.  Were  these  to  be  known  exactly,  or 
exactly  measurable,  the  problem  would  stop  here;  pj  in  (5.256) 
could  provide  the  upper  bound  on  the  error  rate,  arid  p^  in  (5.257) 
the  lower.''  In  the  absence  of  such  complete  statistical  informa- 
tion, one  can  only  develop  error  bounds  by  maximizing  the  integral 
in  (5.256)  and  minimizing  the  integral  in  (5.257).  When  this 
approach  is  taken,  the  incomplete  statistical  information  avail- 
able from  the  measurements  plays  the  important  role  of  constrain- 
ing the  extremal  values  that  either  integral  can  take  on. 

As  in  Section  5 . 7 . 1 , the  first  of  these  constraints  will  be  on 
the  average  value  of  the  interference  variable.  The  second 


Actually,  a single  integral  may  be  sufficient  in  some  situations 
(see  Section  5.7.1),  in  which  case  error  rate  would  be  known 
exactly . 
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will  be  on  p,  the  probability  that  the  interference  variable 
exceeds  some  fixed  threshold  value.  This  differs  from  the  "crest 
factor"  constraint  used  in  Section  5 . 7 . 1 , but  the  work  of  Section 
4.5  indicates  that  it  can  be  more  rapidly  and  reliably  measured. 
The  detailed  forms  of  the  new  set  of  constraints  will  be  presented 
below  in  Sections  5. 7. 3. 3 and  5. 7. 3. 4,  where  we  use  them  to  bound 
probability  of  error  for  nondiversity  operation. 

To  determine  the  upper  limits  of  integration  in  (5.256)  and 
(5.257),  we  make  use  of  the  factor,  7^.  Recall,  from  Section  5.7.1, 
that  this  just  puts  an  upper  limit  on  the  ratio  of  peak  noise 
power  to  average  signal  power  and  arises  because  of  the  dynamic 
limitations  of  receivers;  their  outputs  saturate  when  the  input 
voltage  sufficiently  exceeds  the  average  signal  power  for  which 
the  receiver  was  designed.  Thus,  from  (5.204),  we  have 


Al  - LyL  (UV-  0)  (5.258) 


and 


A 


(5.259) 


In  a similar  vein,  a probability  constraint  on  U given  by 


p = Pr|u>0'UL] 


(5.260) 


where  /?'  is  some  number  between  0 andA  1.  is  exactly  equivalent  to 
the  probability  of  the  event  f2A>/3'Uj  j in  the  mutually  exclusive 
interference  case,  UV=0.  Hence,  forJ  the  integrals  in  (5.256) 
and  (5.257),  the  probability  constraints  will  take  the  form 


p = Pr{A>/S'LyJ  UV  = 0 (5.261) 

and 

p = Pr-[u>  0'2LyJ  (5.262) 

We  close  this  introductory  section  by  demonstrating  that  the 
p constraints  in  (5.261)  and  (5 . 262)  automatically  place  a constraint 
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on  a and,  hence,  the  signs 1-to-noise  ratio.  In  fact,  from  (5.259) 
we  have  immediately 


aL  - U<2LyL 


(5.263) 


Because  of  the  special  constraints  here,  however,  this  is  not  the 
tightest  upper  bound  on  The  tightest  bound  will  be  deter- 

mined below.  First,  however,  we  treat  the  simpler  problem  of 
determining  the  greatest  lower  bound.  Utilizing  the  following 
integral  forms  for  these  constraints 


2LT, 


Q!t 


/ 

0 


UWL(U)  du 


2LyT 

A r L 

p - J WL(U)dU 


and  making  the  change  of  variable 

U 

z ' 2^i 

in  (5.264)  we  have 

o 1 

«L  ■ (2LyL)  J zWL(2LyLz)dz 
0 

, 1 

* (2L vLrfi'f  wL(2LyLz)dz 
p' 


(5.264) 


(5.265) 


(5.266) 


(5.267) 


We  note  that  equality  is  possible  with  an  impulsive  pdf.  Hence 
this  is  the  greatest  lower  bound  on  a^.  By  making  the  same  change 
of  variables  in  (5.265),  the  lower  bound  in  (5.267)  can  be  simpli- 
fied. We  obtain 


aL  & (2LyL)  0'p 
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(5.268) 


Now,  we  find  the  largest  possible  value  of  <*L  subject  to  the 
representations  in  (5.264)  and  (5.265)  and  the  additional  obvious 
constraint 


2Ly  L 

1 - J WL(U)dU  (5.269) 

0 

To  do  this,  we  make  use  of  the  cumulative  distribution,  CL(x), 
where 


x 

CL(x)  » / WL(U)  dU  (5.270) 

0 

We  note  immediately  from  (5.269)  and  (5.265)  that 


CL(2LyL)  - 1 (5.271) 

and 

CL(2/3'LyL)  » 1-p  (5.272) 

Integrating  aL  in  (5.264)  by  parts 

2L7l 

aL  = 2LTl  - J CL(U)dU  (5.273) 

0 


it  becomes  clear  that  we  need  only  minimize  the  integral  in 
(5.273).  Because  of  the  fact  that  C^(x)  may  be  zero  uncil 
x = 2/3/LyL,but  then  must  jump  to  the  value  1-p,  we  note  that 


2LrL 
Min  f 
0 


2L7l 

C^(x)  dx  = Min  J CL(x)  cix 
2j9t rT 

L 


(5.274) 


which,  because  of  the  fact  that  C^(x)  is  monotonically  non- 
decreasing, has  the  value 
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(5.275) 


2L y. 

li 

Min  f CL(x)  dx  - (1  - p)  2LyL(l  - p) 
2j3lyT 

1 J 


Combining  with  (5.273)  we  have  finally  that 


(2LyL)iS/p^  oL  < (2LVL)  [1  - (l-p)(l-^J  (5.276) 

In  terms  of  the  signal-to-noise  ratio,  this  restriction  takes  the 
form 


1 

rLU  - d-p)(i-/J)J  * p 


V'p 


(5.277) 


For  decreasing  /?/and  p,  the  range  of  possible  p values  is  trans- 
lated  to  higher  values  as  one  would  expect. 

5 . 7 . 3 . 3 Error  Rate  Bounds:  Analysis  and  Numerical  Results 

In  this  section  we  develop  the  analysis  leading  to  error 
rate  bounds  for  nondiversity  operation  and  present  the  results  of 
the  numerical  calculations.  Though  only  the  results  for  non- 
diversity operation  are  presented  herein,  we  point  out  that  an 
examination  of  the  earlier  results  (utilizing  different  con- 
straints) and  the  analytical  procedures  indicates  that  the  bounds 
for  diversity  operation  would  be  even  tighter  than  those  pre- 
sented here;  this  is  especially  true  in  the  case  of  the  upper 
bound . 

We  first  deal  with  the  problem  of  determining  the  upper 
bound  on  probability  of  error.  In  keeping  with  the  comments 
accompanying  (5. 256), we  attempt,  through  proper  choice  of  F^(A), 
to  maximize  the  probability  of  error  integral 

n 

pl  » / F1(A)H]_(A)dA  (5.278) 

0 

subject  to  the  constraints 
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(5.279) 


al  m I AF1<A)dA 
0 

and 

yi 

p a j Fj  (A)  dA 


0<£'<1  (5.280) 


and  to  the  additional  obvious  constraint  that  the  candidate  pdf 
integrate  to  1.  Applying  the  completeness  theorem  in  [5.17] 
separately  to  each  of  the  intervals  (0,3/y^)  and  we 

find  that  a general  form  for  the  maximizing  pdf  is 


pUsX(A)  . qifi(A.Ai)  + q26(A-A2)  (5.281) 

where  A^</3*y,  and  A„  ^ /3'V-^  * We  note  that  this  candidate  pdf  has 
the  four  unspecified  parameters,  q^,  A^,  q2,  and  A 2»  Obviously, 

q1  + q2»l  (5.282) 

whereas  application  of  the  constraint  on  in  (5.279)  gives 


“l  = qlAl  + q2A2  (5.283) 

and  application  of  the  p constraint  in  (5 . 280)  just  indicates  that 
q2ap,  so  that  now  the  candidate  pdf  has  just  one  unspecified 
parameter,  A^;  i.e., 

F^aX(A)  - (l-p)6(A-A1)  +p6(A-A2)  (5.284) 

where 

O'.  - (1  - p)A1 

A = — (5.285) 

2 p 
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Thus,  the  probability  of  error 


Max 

P1 


(l-p)H1(A1)  + pH1(A2) 


(5.286) 


must  be  maximized  with  respect  to  the  unspecified  parameter  A^. 
To  accomplish  this,  we  note  that 


dp 


Max 


d A. 


|a  - p) 


(5.287) 


and  that  the  second  derivative  is  negative  throughout  the  inter- 
val (a  consequence  of  the  fact  that  A2>A^).  It  is  immediately 
apparent  from  (5.287)  that  p^ax  is  maximized  at  A-^  = A2«  Thus, 
from  (5  c 285)  we  determine  that  the  maximum  occurs  at  A^=a-^;  the 
upper  bound  on  probability  of  error  has  the  value 


Max  _ „ , x 

pi  - W 


(5.288) 


This  is  the  same  as  the  upper  bound  determined  under  the  earlier 
constraint  problem  [see  (5.220)]. 

We  now  turn  to  the  problem  of  determining  the  lower  bound  on 
the  probability  of  error.  In  keeping  with  the  comments  accom- 
panying (5.257),  we  attempt  through  proper  choice  of  W]i(A),  to 
minimize  the  probability  of  error  integral 

2yl 

pl  « J W1(U)G1(U)dU  (5.289) 

0 

subject  to  the  constraints 

UWL(U)dU  (5.290) 


Wx(U)dU  (5.291) 
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2 y. 


al  = / 


and 


0 


7 


J 


and,  as  before,  the  obvious  restriction  that  the  pdf  integrate 
to  1.  Applying  the  completeness  theorem,  we  have  as  a general 
form  for  the  minimizing  pdf 


W^in(U)  - q^OJ)  + q2§  (U  - 2£/V1)  + q3fi(U-U1>  (5.292) 

which  has  the  four  unspecified  parameters,  q.  , q_,  q„,  and  U1 . 
Obviously, 

qjL  + q2  + q3  " 1 (5.293) 

whereas  application  of  the  constraint  on  in  (5.290)  gives 

al  = 92(2*/yl)+  q3U  (5.294) 

and  application  of  the  probability  constraint  in  (5.291)  gives 
the  result 

p ■ q2  + q3  (5.295) 

Solving  these  three  equations,  we  obtain  the  results 


qx  = 1 - p 

(5.296) 

p°r  “i 
q2  U1  - 2p'y1 

(5.297) 

a - ( 2f$'y . ) p 

„ J.  1 

q3  U1  - 2/J'y. 

(5.298) 

This  procedure  leaves  one  degree  of  freedom  through  dependence  on 
the  variable  U^.  We  now  minimize 
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.A 


min 

P1 


G^)  (5.299) 


G1(28/y1)  + 


<*!“  (2 P\)p 

lIl“2j8/y1 


with  respect  to  U.  . This  is  best  accomplished  by  first  rearrang- 
ing the  terms  in  (5.299).  We  obtain 


min 

pl 


“ P 


- (28/y1)G1(u1) 


ul  ’ 


+ a. 


G1(U1)  "Gl(2^'yi) 


Ux  - 2j8'y1 


(5.300) 


Since  GX(U)  is  everywhere  concave  downward,  a simple  graphi- 
cal construction  indicates  clearly  that  the  second  bracketed  term 
in  (5.300)  is  minimized  at  Ik  = 2y...  Evaluating  the  derivative  of 
the  first  term,  we  obtain 


d[-J. 

dlL 


2/3  'y. 


(u1-2^'y1) 


G(U  ) -G(2 py  ) 

G,(ui>  +-^-2^y1 


(5.301) 


which,  because  U-,  is  greater  than  2/8  y.  , indicates  that  the  first 
term  in  (5.300)  is  minimized  at  U1=2y  i,  and  q q and  q„  in 

(5.296)  through  (5.298)  take  on  the  values  12  5 


qx  - i - p 

_ (27x)p  - 
q2  2yx(l  - p) 

_ a1  - (2/3,yi)  p 
q3  2y1(i  - /3) 

Turning  now  to  the  dependence  on  a.,  we  note  from  (5.276) 
that  is  restricted  because  of  the  constraints  to  values 
satisfying 


(5.302) 

(5.303) 

(5.304) 
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(5.305) 


al  * ( 2yl ^ p 


Hence,  there  is  no  value  of  that  allows  q^  in  (5.298)  to  take 
on  negative  values.  The  situation  is  not  as  simple  in  the  case 
of  the  second  impulse  which  vanishes  as  reaches  and  exceeds 
(27-.  )p.  In  this  case,  the  only  way  that  the  constraints  can  be 
satisfied,  especially  the  requirement  for  larger  , is  that  the 
lower  impulse  begins  to  move  to  the  right.  The  minimizing  pdf  now 
has  the  form 

W^OJ)  » (l-p)6(U-U,)  + p6(U-2y1)  (5.306) 


where  the  impulse  strengths  follow  directly  from  the  probability 
constraint  and  the  requirement  that  the  pdf  integrate  to  1. 
Applying  the  constraint  on  a^,  it  is  easy  to  see  that 

a,  - p(2 71  ) 

U*  = — (5.307) 

1 - p 

and,  hence,  that 

min 
P1 


al>2yip  (5.308) 

It  would  seem  at  first  glance  that  the  behavior  of  (5.308) 
at  values  of  U1  in  excess  of  2/37  would  require  separate  considera 
tion;  this  is  because  both  impulses  would  fall  in  the  region 
(2/3y^,2 7-)  over  which  the  integral  defining  the  probability  con- 
straint takes  place.  We  note,  however,  that  a value  of  U1  ex- 
ceeding 2/37  implies,  from  (5.307),  a value  of  a ^ satisfying 


= (1  -p)  G, 


- p(271)’ 
1-  p 


+ pG1(271) 


<x1  > P ( 27 x ) + 2j3tl  -P)71 

which  can  be  rearranged  to  take  the  form 


(5.309) 
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(5.310) 


> 2r1[l  - (l-p)d-rt] 


This  places  a,  outside  of  the  range  of  possible  values  given 
in  (5.276).  Thus,  the  lower  bound  on  probability  of  error  has 
been  determined  for  all  possible  values  of  a.,  . 

Recalling  from  (5.202)  that 


p “ 2 /a1 


(5.311) 


we  summarize  the  above  results  in  terms  of  the  signal-to-noise 
ratio: 


Max 


min 


• 


y,  \ 


[\-  2P-. 


+ PG1(2r1) 


1 1 

^[l  - (i-p)(l-^r  < p < rlP 


(5.312) 


(5.313) 


min 


V-p 


y\(i  - ¥) 


G1(2j3'y1)  + 


- - B'r,  p 

9. LI 

y^i  - P) 


Gi<2V 


——  < p < —A — 

rip  pyi p 


(5.314) 


Probability  of  error  bounding  curves  determined  from  numeri- 
cal calculation  of  these  expressions  are  shown  in  Figures  5.34- 
5.37.  The  curves  shown  hold  exactly  for  the  binary  FSK  case  but 
from  the  discussion  in  Section  5 . 7 . 1 can  be  thought  of  as  worst- 
case  curve  both  in  terms  of  probability  of  error  and  tightness 
of  bounds . 

The  first  three  curves  are  for  increasing  values  of  p, 
p = 0.001,  p = 0.01,  and  p = 0.5,  and  various  values  of  threshold  j3'. 
The  increasing  values  of  p imply  increasing  values  of  percentage 
measurement  error  of  p,  as  discussed  in  Section  5. 7. 3.1.  These 
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have,  respectively,  the  approximate  values  ,/1000/M,  JlGOjK,  and 
v/T7m  ( exact),  where  M is  the  number  of  independent  samples,  as 
discussed  in  Section  4. 5. 4. 3. 

In  each  of  these  curves,  the  breakpoint  [see  (5.313)  and 
(5.314)]  occurs  at  l/y,p  which, is  too  far  to  the  left  to  be 
illustrated  in  Figure  5.36. 

In  Figure  5.37,  the  curves  are  plotted  for  fixed  and  various 
F • 

These  curves  indicate  clearly  that  measurements  of  the  proba- 
bility with  which  the  noise  variables  exceed  some  threshold  can 
be  used  to  develop  bounding  curves  for  probability  of  error.  The 
measurements  are  easy  to  obtain  and  with  suitable  operational 
choices  can  be  obtained  with  high  accuracy. 

With  a greater  number  of  thresholds,  the  measurements  would 
correspond  to  a simple  histogram,  and  there  is  no  doubt  that  one 
could  develop  considerably  tighter  bounds. 

5.7.4  Summary 

In  this  section  we  discussed  the  way  in  which  the  inter- 
ference measurements,  covered  earlier  in  Section  4.5,  can  be  used 
to  develop  bounds  on  the  probability  of  error.  To  accomplish  this, 
we  drew  heavily  on  results  published  earlier  in  [5.17].  The  way  in 
in  which  the  error  rate  and  the  error  rate  bounds  are  related  to 
the  measurements  was  discussed  in  Section  5.7.1,  In  this  section, 
we  demonstrated  how  the  probability  of  error  integrals  can  be  ex- 
tremized  by  applying  constraints  on  average  interference  power  and 
on  the  interference  crest  factor.  The  fact  that  these  constraints 
can  lead  to  tight  bounds  on  error  probability  was  evidenced  by  the 
curves  presented  in  Figures  5.29-5.31  which  were  taken  from  [5.17]. 
In  Section  5.7.2  we  discussed  the  accuracy  with  which  the  error 
rate  bounds  can  be  determined  from  the  interference  measurements; 
as  one  would  expect,  the  errors  in  developing  the  bounding  curves 
depend  in  direct  fashion  on  the  errors  in  the  measurement  of  the 
interference  mean  value,  maximum  value,  and  on  the  error  in  esti- 
mating the  mean-squared  value  of  the  channel  transfer  function. 

The  relationship  between  the  rins  errors  in  estimating  these  quan- 
tities and  in  estimating  the  bound  values  was  given  in  Section 
5. 7. 2. 2 [see  (5.254)]  and  illustrated  in  Figure  5.33.  It  is  clear 
from  the  work  of  this  section  that  reliable  interference  measure- 
ments can  lead  to  reliable  bounds  on  error  probability. 

In  Section  5.7.3  we  tried  an  alternate  approach  to  the  prob- 
lem of  developing  bounds  on  probability  of  error.  This  approach 
made  use  of  estimates  of  p,  the  probability  that  the  interference 
variable  U exceeds  some  fixed  threshold  U1  (see  Section  4. 5. 4. 4). 
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Thus,  the  measurements,  in  effect,  make  up  a simple  two  bin  histo- 
gram for  the  interference.  It  seems  reasonable  to  expect  that 
variants  of  this  approach  could  be  used  to  bound  error  probability 
when  histograms  with  a higher  number  of  bins  are  available. 
Although  the  accuracy  of  determining  the  bounds  was  not  investi- 
gated as  thoroughly  as  that  determined  earlier,  the  indication  in 
Section  4. 5. 4. 4 that  p can  be  measured  reliably  indicates  also 
that  the  bounds  can  be  estimated  reliably. 
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SECTION  6 

TECHNIQUES  UTILIZING  SPECIAL  PROBING  SIGNALS 


The  previous  sections  have  been  concerned  with  the  area  of 
main  interest  to  the  present  study,  the  in-service  estimation  of 
modem  performance  based  upon  propagation  medium  limitations  alone 
(i.e.,  assuming  a well-functioning  modem)  without  the  use  of 
special  probing  signals.  A number  of  concepts  were  developed 
and  analyzed  yielding  a number  of  techniques  of  wide  applicability. 
However,  there  are  cases  where  special  probing  signals  are  needed 
for  the  most  reliable  performance  estimation.  Thus,  in  the  case 
of  LOS  channels  where  deep  fades  can  sometimes  be  highly  frequency- 
selective,  the  complex  transfer  function  T(f,t)  is  needed  to  esti- 
mate the  performance  of  conventional  suboptimum  modems  which  do 
not  attempt  maximum  likelihood  demodulation. 

It  is  not  possible  to  estimate  a fading  dispersive  channel's 
delay  power  spectrum  Q (4 ) from  |T(f,t)|2.  The  most  accurate  pre- 
diction of  irreducible  error  rate  for  time-gated  and  anti-multipath 
troposcatter  modems  is  with  the  use  of  Q(£).  The  measurement  of 
Q(£)  requires  the  use  of  special  probing  signals  if  processing  is 
to  take  place  at  RF  or  IF. 

The  HF  channel  was  found  to  be  the  most  difficult  to  esti- 
mate modem  performance  using  the  techniques  of  the  previous 
sections.  Special  probing  signals  could  reduce  estimation  error. 

The  signal  distortion  properties  of  the  channel  are  sometimes 
characterized  by  measurements  at  three  levels  of  decreasing 
knowledge: 

1)  Measurement  of  system  functions  (e.g.,  time-variant 
transfer  function  and  impulse  response). 

2)  Measurement  of  correlation  functions  of  system 
functions  (e.g.,  frequency  correlation  function, 
delay  power  spectrum,  etc.). 

3)  Measurement  of  gross  parameters  of  system  function 
(e.g.,  coherence  bandwidth,  Doppler  spread,  etc.). 

We  consider  each  of  these  in  order. 
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6.1  Out-of-Service  Measurement 


Here  we  assume  that  service  has  been  interrupted  and  it  is 
desired  to  measure  the  channel  characteristics. 

6.1.1  System  Function  Measurement 

The  problem  of  probing  and  measuring  system  functions  for 
time-variant  channels  has  been  studied  extensively  by  Bello  and 
Esposito  [6.1],  Bello  [6. 2], [6. 3],  and  Bello  and  Pinto  [6.4]. 

It  has  been  found  that  for  well-designed  signal  processing, 
the  accuracy  in  measurement  of  the  system  functions  is  essentially 
invariant  to  the  technique  used.  The  system  function  channel 
measurement  technique  used  most  often  is  the  cross-correlation 
technique  employing  a periodic  probing  signal  containing  a pseudo- 
noise sequence  of  + 180° modulating  a carrier.  A simplified  block 
diagram  of  the  receiver  processing  is  shown  in  Figure  6.1.  The 
received  multipath-distorted  PN  waveform  is  mixed  with  a set  of 
delayed  replicas  of  the  transmitted  sequence  to  extract  samples 
of  the  channel's  complex  time-variant  impulse  response  at  the 
outputs  of  the  lowpass  filters.  The  filters  have  passbands  flat 
over  the  spectrum  occupied  by  the  fading  channel,  i.e.,  the 
"total"  Doppler  spread. 

If  z(t)  denotes  the  complex  envelope  of  the  transmitted 
probing  signal,  g(t,£)  the  complex  envelope  of  the  impulse  re- 
sponse of  the  channel,  and  rj(t)  the  complex  envelope  of  the 
received  additive  noise,  then  the  received  probing  signal  is 
given  by 


w(t)  33  J*  z(t  - £)  g(t,£)  d£  + T7(t)  (6.1) 

The  cross-correlation  operation  to  produce  a (time- varying) 
sample  of  the  impulse  response  at  £~t,  g(t,r),  is  given  by 

g(t,r)  - [w(t)z*(t  - r )]  ® k(t)  (6.2) 

where  k(t)  is  the  impulse  response  of  the  lowpass  filter.  For 
this  discussion,  four  assumptions  will  be  made: 
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Figure  6.1  Simplified  Block  Diagram  for  Correlation  Technique  of  Impulse  Response 


1)  The  probing  signal  power  spectrum  is  flat  over  W, 
the  bandwidth  of  the  channel  that  it  is  desired  to 
characterize. 

2)  The  lowpass  filter  k(t)  is  nondistorting  to  the  time- 

variant  impulse  response  sample  fluctuation.  This 

means  that  the  nondistorting  bandwidth  of  k(t)  is 

equal  to  or  greater  than  the  total  Doppler  spread 

of  the  channel,  EL  ^ • 

tot 

3)  The  bandwidth  W is  much  bigger  than  B . 

4)  The  product  B L»  1>  where  L is  the  duration  of  the 
channel  impulse  response. 

Assumptions  3)  and  4)  are  valid  for  the  channels  under  dis- 
cussion. Assumptions  1)  and  2)  are  convenient  for  analysis 
although  not  optimum  for  minimizing  the  estimation  error  due  to 
the  combined  effects  of  additive  noise  and  distortion.  However, 
it  will  be  seen  below  that  by  comparison  with  the  incoherent 
probing  techniques  of  the  previous  sections,  the  output  additive 
noise  is  very  much  smaller  so  that  we  may  be  more  generous  here 
in  allowing  some  extra  noise  through,  to  essentially  remove  dis- 
tortion effects. 

With  the  above  assumptions,  the  estimated  impulse  response 
becomes 


g(t ,r)  “ J*  C (r  - O g(t,£)d£  + 77 (t ,t ) (6.3) 

£* 

where  G (r)  is  the  autocorrelation  function  of  the  probing  signal, 
and  r)(t7,T ) is  a noise  (in  effect  a ’’noise  channel")  with  power 
spectral  shape  (for  fixed  r)  proportional  to  |K(f)|2. 

All  radios  have  filters  at  the  transmitter  and  receiver 
which  limit  bandwidth.  These  may  be  defined  to  be  part  of  the 
channel  to  be  probed  or  not,  depending  upon  the  objectives  of 
the  probing.  In  either  case,  the  probing  signal  would  be  de- 
signed to  have  a flat  spectrum  over  the  band  of  interest  and 

h(t ,t)  - J C (t  -O  g<t,$)d?  (6.4) 

z 

would  be  the  desired  impulse  response. 


It  may  be  shown  from  the  referenced  work  by  Bello  that,  for 
a dispersive  channel  with  delay  power  spectrum  Q(£),  the  strength 
of  the  desired  term  to  the  estimation  mean-squared  error  is 
is 


Pout  W Pin  B 


K 


;Q(4)d£ 


(6.5) 


where  p,  is  the  input  SNR  measured  in  the  bandwidth  of  the 
received  signal.  For  a uniform  delay  power  spectrum  of  width  L, 
(6.5)  becomes 


pout  ^ pin  B^L 
K. 


(6.6) 


Since  Bj^  can  be  set  close  to  Btot>  the  total  Doppler  spread  of  the 
channel,  and  since 


s s B.  . L «<  1 
tot 


(6.7) 


for  the  channels  of  interest,  we  see  that  excellent  measurement 
accuracy  can  be  obtained  with  out-of- service  probing  because  the 
communication  system  is  normally  designed  so  that  p^n  is  » 1. 

6.1.2  Correlation  Function  and  Gross  Parameter  Measurement 

Since,  on  an  out-of-service  basis,  the  system  function  can 
be  measured  with  high  accuracy  for  the  channels  of  interest,  the 
only  significant  contribution  to  the  estimation  error  for  measure- 
ment of  a channel  correlation  function  will  be  the  channel  fluc- 
tuations which  must  be  averaged  out.  One  may  show  that  the  rms 
fractional  estimation  error  a always  takes  the  form 

6 « (6.8) 

JK 

where  K is  an  effective  number  of  independent  samples  of  the 
channel  during  the  averaging  time.  This  may  be  equated  approxi- 
mately to  the  product  of  the  observation  time  and  the  rms  Doppler 
spread  of  the  channel.  Thus,  for  example,  if  it  is  desired  to 
measure  the  delay  power  spectrum  of  a channel  Q(£),  the  rms 
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estimation  error  due  to  channel  fluctuations  can  be  expected  to 
be 


(6-9) 

a/K 

Of  course,  if  it  is  desired  to  measure  Q (4 ) where  it  is  very 
small,  errors  due  to  the  additive  noise  would  have  to  be  con- 
sidered. From  (6.5)  - (6.7)  we  see  that  as  long  as 


» — (6.10) 

Tnax  pin 


the  noise  effects  may  be  disregarded. 


The  same  three  classes  of  techniques  may  be  used  for  gross 
channel  parameter  measurement  with  special  probing  signals  a3 
were  discussed  in  Section  4 for  the  received  signal  probe  case. 

In  the  present  case,  coherent  processing  techniques  may  be  used 
for  the  differentiation  and  level  crossing  approaches,  and  an 
improved  correlation  function  estimate  may  be  used  in  the  correla- 
tion technique.  Moreover,  in  addition  to  rms  Doppler  spread,  a 
mean  Doppler  shift  may  be  measured.  We  briefly  consider  these 
measurement  techniques. 


If  a carrier  is  transmitted,  the  received  process  v/hen 
referenced  to  the  carrier  frequency  has  the  spectrum  P(f).  It 
is  convenient  to  deal  with  the  complex  envelope  of  this  received 
process  z(t),  which  has  just  this  low-frequency  spectrum  P(f). 

In  [6.1]  it  is  shown  that  the  Doppler  spread  parameter  B is  given 
by  the  following  operations  on  z(t): 


B = 


IT 


| < m2> 

<z*z> 

||  < !z(t)|  2> 

<lz(t)|2> 

(6.11) 


where  the  triangular  brackets  indicate  time  averages. 

In  actuality,  the  operations  upon  z in  (6.11)  would  most 
likely  be  carried  out  by  operating  on  the  "in-phase"  and 
"quadrature"  components  x(t)  and  y(t)  related  to  z(t)  by 
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z(t)  - x(t)  + jy(t ) 


(6.12) 


It  is  shown  by  Bello  [6.5)  that  in  terms  of  x(t)  and  y(t), 
B is  given  by 


B 


.Ui>2  + ,(±)_2>. 

(xy- yx) 

/ <*2  + y2> 

<x  + y ) 

(6.13) 


The  in-phase  and  quadrature  components  (x,y)  can  be  deter- 
mined by  multiplying  the  received  carrier  by  both  a local  carrier 
and  a 90°  shifted  local  carrier  at  the  same  frequency  as  the 
received  carrier  (or  as  near  to  the  same  frequency  as  possible) 
and  then  extracting  the  low-frequency  components.  Strictly 
speaking,  D is  independent  of  mean  Doppler  shift,  and  thus  pre- 
cise knowledge  of  the  received  carrier  frequency  is  not  necessary. 
However,  as  the  local  carrier  frequency  departs  from  the  received 
carrier  frequency,  the  extracted  x(t)  and  y(t)  increase  in  band- 
width, necessitating  larger  bandwidth  filters  and  passing  more 
noise.  Thus,  from  the  point  of  view  of  maximizing  signal- to- 
noise  ratio,  it  is  desirable  to  keep  the  local  carrier  frequency 
as  near  as  possible  to  the  received  signal  frequency.  An  esti- 
mate of  the  mean  Doppler  shift  defined  as  the  centroid  of  the 
power  spectrum  is  provided  by 


f 


_L 

2n 


( xy  - yx) 

2 2 
(x  + y ) 


(6.14) 


In  practice,  the  derivatives  may  be  approximated  by  differences 
as  was  done  for  the  envelope  techniques  in  Section  4. 

The  measurement  of  multipath  spread  by  differentiation  tech- 
niques follows  an  approach  which  is,  in  principle,  dual  to  that 
described  above  for  Doppler  spread,  where  the  duality  involves  an 
interchange  of  the  roles  of  time  and  frequency.  The  frequency 
derivatives  would  be  approximated  by  differences  between  the 
complex  amplitudes  of  two  received  tones  judiciously  spaced  in 
frequency. 

For  the  level  crossing  approach,  we  may  show  from  the  results 
of  Rice  [6.6]  that  the  average  number  of  times  m the  in-phase  or 
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quadrature  components  of  narrowband  C.uassian  noise  crosses  zero 
is  given  identically  by  the  rms  Doppler  spread,  i.e., 

m “ B (6.15) 

Thus,  by  counting  zero  crossings  in  some  finite  time  interval,  we 
obtain  an  estimate  of  the  rms  Doppler  spread. 

Also,  in  an  exactly  analogous  fashion,  m,  the  number  of 
times/hertz  that  the  real  (or  imaginary)  part  of  the  channel 
transfer  crosses  zero  is  identically  equal  to  the  rms  multipath 
spread  L,  i.e., 

m “ L (6.16) 


The  correlation  techniques  for  measurement  of  multipath  and 
Doppler  spread  were  discussed  in  Section  4.  In  that  section, 
they  were  based  upon  the  channel  envelope-squared  time  or  fre- 
quency correlation  functions.  Here  the  same  procedure  would  be 
used  but  they  would  be  based  upon  the  correlation  function  of  the 
real  or  imaginary  parts  of  the  transfer  function. 

The  estimation  errors  for  multipath  and  Doppler  spread  would 
also  be  dependent  primarily  on  the  number  of  independent  channel 
samples  in  the  averaging  interval.  The  rms  fractional  error  6 
would  vary  as  1/VK~,  i.e.. 


6 ^ 


(6.17) 


6 . 2 In-Service  Measurement  Techniques 

The  previous  section  considered  various  classes  of  channel 
measurement  techniques  utilized  in  the  absence  of  data  signal 
transmissions.  Out-of-service  operation  is  highly  undesirable. 
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One  may  obtain  channel  measurement  information  utilizing  probing 
signals  in  the  same  frequency  band  and  during  in-service  opera- 
tion by  multiplexing  the  data  signal  and  probing  signal  together. 
We  may  consider  four  basic  multiplexing  arrangements: 

1)  Data-probing  time-division  multiplex 

2)  Data-probing  frequency-division  multiplex 

3)  Data-probing  FDM/TDM 

4)  Nonorthogonal  multiplex 

In  each  of  these  arrangements,  the  objective  is  to  multiplex 
and  demultiplex  the  data  and  probing  signals  on  an  essentially 
noninterfering  basis.  In  the  case  of  time-division-multiplex, 
probing  signals  are  inserted  in  series  with  data  signals.  Note 
that  the  method  of  multiplexing  incoming  digital  data  can  gener- 
ally be  independent  of  the  way  the  probing  and  data  signals  are 
multiplexed  together.  Thus,  the  TDM  (time-division-multiplex) 
of  probing  and  data  signals  can  occur  with  either  FDM  or  TDM  of 
data  subchannels.  Furthermore,  the  multiplexing  of  probing  and 
digital  data  can  occur  at  baseband,  RF,  or  IF,  yielding  a large 
number  of  potential  combinations  to  be  examined. 

The  concepts  of  TDM,  FDM,  and  TDM/FDM  are  best  understood  by 
imagining  that  the  available  time-frequency  space  (either  at 
baseband,  IF,  or  RF)  is  divided  up  into  two  mutually  exclusive 
regions  assigned  to  the  probing  signal  and  to  the  data  signal. 

In  order  for  there  to  be  little  loss  in  channel  capacity,  the 
percentage  of  time-frequency  space  devoted  to  probing  signals 
must  be  a small  fraction  of  the  total  available  time -frequency 
space . 

The  first  three  techniques  are  orthogonal  multiplexing 
arrangements.  The  last  technique  has  the  advantage  of  using  no 
additional  time-frequency  space.  However,  it  cannot  always  be 
used  because  of  the  mutual  interference  created.  We  explore  this 
point  for  the  use  of  an  in-band  wideband  probing  signal.  There 
are  two  requirements  for  utilization  of  this  technique: 

1)  The  probing  signal  must  be  at  a low  level  relative 
to  the  data  signal  at  transmission,  typically  20  dB 
or  lower  to  avoid  degrading  the  communications. 

2)  The  processing  gain  of  the  probing  signal  demodulator 
must  be  high  enough  to  overcome  the  large  interfering 
data  signal. 
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Generally  speaking,  the  probing  technique  can  be  any  wideband 
coherent  processing  technique  which  spreads  the  probing  signal 
over  time- frequency  space.  However,  care  must  be  exercised  in 
selection  of  this  technique  to  see  that  discrete  spectral  com- 
ponents are  either  low  enough  or  properly  placed  so  as  not  to 
produce  a local  region  of  low  data/prober  power  density  ratio. 
Presuming  this  care  has  been  taken,  we  can  estimate  from  Eqs. 
(6.5)-  (6.7)  what  the  relation  between  channel  parameters  must  be 
to  allow  use  of  the  nonorthogonal  multiplexing  technique  for 
measuring  channel  system  functions.  In  this  equation,  we  inter- 
pret p^n  as  the  ratio  of  prober  power  to  data  power,  a value  we 
have  quoted  as  being  desirably  -20  dB  or  smaller.  If  we  desire 
at  least  a 20-dB  measurement  SNR,  pout,  then  from  Eqs.  (6.6)  and 
(6.7)  we  see  that  the  spread  factor 


s £ 10“4  (6.18) 

For  troposcatter  links,  multipath  spreads  are  not  much  bigger 
than  1 /is,  while  Doppler  spreads  range  from  a fraction  of  1 Hz 
to  10  Hz.  Thus,  (6.18)  will  be  well-satisfied  for  troposcatter 
links.  In  the  case  of  microwave  LOS  links,  s is  even  smaller, 
rarely  exceeding  10~10  to  io-ll.  Thus,  LOS  links  are  ideally 
suited  for  in-service  probing.  A similar  calculation  for  iono- 
spheric scintillation  channels  reveals  these  suitable  also.  For 
HF  channels,  multipath  spreads  are  usually  less  than  2 ms  and 
Doppler  spreads  less  than  1 Hz.  This  yields  s = 2 x 10-3  which 
violates  (6.18).  Thus,  some  difficulty  may  be  anticipated  in 
utilizing  a low-level  probing  signal  in  parallel  with  data  trans- 
mission for  in-service  probing  on  HF  links.  In  the  case  of  HF 
links  which  are  SSB,  however,  one  or  more  frequency  slots  are 
usually  provided  for  pilot  tones.  These  may  be  used  for  accurately 
measuring  gross  channel  parameters  and  the  complex  transfer  func- 
tion at  those  frequencies. 

From  a complexity  point  of  view,  it  would  be  best  to  effect 
in-service  probing  in  a way  that  would  involve  least  changes  in 
existing  equipment.  The  two  techniques  that  make  most  sense  from 
this  point  of  view,  when  applicable,  are  nonorthogonal  multiplex 
with  low-level  wideband  probing  signals  and  tones  placed  at  the 
edge  of  or  in  the  signal  band.  Of  course,  a series  of  tones 
placed  across  the  band  may  be  regarded  as  a wideband  probing 
signal.  We  consider  the  applicability  of  these  techniques  for 
each  channel. 
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In  the  case  of  HF  channels,  we  have  already  pointed  out  that 
a low-level  PN  prober  would  not  be  satisfactory  because  the  spread 
factor  of  the  channel  is  too  large.  Since  SSB  with  parallel  sub- 
carriers is  used  in  HF,  it  is  a simple  matter  to  insert  probing 
tones  in  place  of  data  subcarriers.  This  would  always  be  at  the 
cost  of  reducing  capacity,  but  frequently  at  least  one  tone  is 
provided  anyway  for  AFC. 

For  the  other  channels,  transmission  signals  are  constant 
envelope  and  nonlinear  essentially  hardlimiting  power  amplifiers 
are  used.  In  this  connection  we  note  that  if  a low-level  probing 
signal  p(t)  plus  a constant  envelope  data  signal  exp[j<p(t)]  are 
fed  to  a hard  limiter,  the  output  may  be  expressed  in  the  form 


<>(t)  + n(t) 

| ) 4.  p(t)  | 


= ejo(t)  + |p(t)  - “p*(t)ej2<p(t)  + 


(6.19) 


2 

where  the  missing  terms  are  of  the  order  of  p . By  keeping  the 
strength  of  p(t.)  20  - 30  dB  below  the  constant  envelope  data  signal, 
tile  additional  terms  may  usually  be  neglected.  The  only  proviso 
is  the  possible  creation  of  discrete  frequency  components  which, 
although  small,  may  fall  on  critical  parts  of  the  data  signal 
structure.  This  possibility  must  be  examined  in  each  case. 

.suming  that  the  other  terms  may  be  neglected,  we  see  that 
the  „ ^obing  signal  has  been  reduced  6 dB  by  the  hard  limiting  and 
a spurious  term  comparable  in  strength  to  the  probing  signal  has 
appeared.  The  spurious  term  will  only  be  harmful  to  the  channel 
measurement  if  expf  j2<a(t)  j is  not  a broadband  signal.  To  take  an 
extreme  case  if  the  data  signal  is  binary  PSK,  exp[j2<p(t)]  would 
be  a constant  and  the  spurious  signal  would  be  the  negative  con- 
jugate of  the  probing  signal.  Barring  this  special  case,  the 
spurious  signal  would  be  a wideband  signal  contributing  much 
smaller  data  noise  to  the  channel  measurement  than  the  data  sig- 
nal itself.  In  all  cases  of  interest,  the  spurious  signal  pro- 
duces no  harm  to  the  data  signal  demodulation  because  it  is  of 
the  same  strength  as  the  output  probing  signal  which  is  already 
at  a very  low  level. 

In  view  of  the  above,  it  would  appear  that  either  PN  or 
parallel  in-band  probing  is  feasible  for  LOS  and  troposcatter 
links,  the  ultimate  choice  being  cost. 
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For  satellite  ionospheric  scintillation  channels,  multipath 
is  not  a factor  of  importance  and  a single  low-level  tone  placed 
adjacent  to  the  data  signal  can  be  used  to  characterize  the 
channel. 

6.3  Error  Rate  Estimation 


The  error  rate  estimation  techniques  discussed  in  Section  5 
based  upon  processing  the  estimated  values  of  the  squared  magnitude 
of  the  transfer  function  may  obviously  be  used  when  the  actual 
complex  transfer  function  or  impulse  response  is  known.  Moreover, 
the  bias  errors,  which  were  especially  troublesome  for  HF  channels, 
will  vanish  and  the  variances  of  the  estimation  errors  will  reduce. 
The  latter  reduction  will  not  be  so  dramatic  as  the  reduction  of 
biases  because  the  unaveraged  channel  fluctuations  usually  repre- 
sent the  major  source  of  estimation  error  variance. 

Aside  from  the  general  benefits  for  HF  channels,  perhaps  the 
biggest  bonus  comes  from  the  improved  ability  to  estimate  the 
irreducible  error  probability  produced  by  excessive  multipath  in 
high-speed  troposcatter  modems  and  the  condition  error  probability 
due  to  deep  frequency-selective  fades  on  LOS  links.  Insufficient 
time  was  available  to  evaluate  the  effectiveness  of  such  estima- 
tion techniques. 
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APPENDIX  A 


GROSS  CHANNEL  PARAMETER  ESTIMATION  FROM 
BIAS  UNCORRECTED  ESTIMATES 


A.l  Introduction 


In  this  appendix  we  address  the  problem  of  estimating  gross 
channel  parameters  from  bias  uncorrected  estimates  of  the  magni- 
tude squared  of  the  channel  transfer  function.  In  particular, 
we  will  estimate  the  rms  Doppler  and  multipath  spreads  using  the 
differentiation  technique  considered  in  Section  4.3,  and  the 
mean  SNR. 

It  will  be  shown  that  using  bias  uncorrected  estimates  does 
not  appreciably  affect  the  performance  of  these  gross  channel 
parameter  estimators.  Thus,  the  above  gross  channel  parameters 
can  be  accurately  estimated  without  bias  correction,  and  then 
used  to  correct  the  bias  in  the  magnitude-squared  channel  trans- 
fer function  estimates. 


When  the  bias  of  the  magnitude  squared  of  the  channel  trans- 
fer function  estimate  is  not  removed,  this  estimate  can  be 
expressed  as  (see  Section  4.1) 


|T(F,tp)|2  - |T(F,tp)|2(l  + ep)  + 6p 


where 


E{ep}  = -a  = - HL2  - KB2 

E{6p}  - a E{|T(f,t)|2}  + 2N0/Pz(F) 


CgE2{|T(f,t)|2} 

2 


(1) 


(2) 


2 

In  the  above  expressions,  ]T(F,t  )|  is  the  estimate  of  the 
magnitude  squared  of  the  channel  tranlfer  function , j T (F, t )|^  is 
its  actual  value,  and  € and  are  given,  respectively, 

(4.57)  and  (4.100).  b 
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A. 2 Estimation  of  the  Mean  SNR 

The  estimation  of  the  mean  SNR  is  performed  by  using  the 
estimates  of  the  magnitude  squared  of  the  channel  transfer  func- 
tion to  form  estimates  of  the  instantaneous  SNR,  which  are  then 
averaged  to  give  the  mean  SNR  estimate.  We  can  define  an  instan- 
taneous SNR  by 


P (F)  , 

*V  “is7  iT<F.yr  <3) 

and  a mean  SNR  by 

r - E{y(tp)}  (4) 

Using  (1)  we  can  estimate  T by 

- P (F)  t KT  9 1 

r - isj-  I gL'^’V1  <1+V  + 6PJ  <5> 


Averaging  over  the  errors  and  the  channel  fluctuations  gives 


r + i 


(6) 


where  the  estimation  errors,  c and  6 , were  assumed  to  be  inde- 
pendent of  the  channel  fluctuations.  ^ 


Therefore,  no  correcting  for  the  bias  in  estimating 
|T(F,t  )|2  produces  a bias  in  estimating  T.  However,  since  this 
bias  is  a known  constant,  it  can  be  corrected. 


The  second  moment  of  the  mean  SNR  estimate  is  given  by 
o2 


E{f2} 


J 


PZ(F) 


|L  2NoJ 


? s U' 


T(F,tp)|  (l  + «p)  + 6p 


iT(F,tp)ra+fq)  + sq 


(7) 
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After  some  simplification,  we  can  express  E-lT  \ by  (where  a 
complex  Gaussian  scatter  channel  is  assumed) 


f2}  * \ 


2T 


(1  -a)  + e +• 


11 : a)2!'2  v 


K 


K2 


S.  q?L  (1+PP<1) 


+ (oT-tl)z  4-  2T(1  - a)  (aT+  1) 

2 


(8) 


p=l  q 

pj^q 


where  p is  the  magnitude  of  the  normalized  channel  correlation 
function^and  is  defined  by 


A |E{T*(F,tp)T(F,tq)} 
E{|T(F,t)!2} 


(9) 


A 

* Using  the  mean  of  T as  given  byA(6)  and  the  second  moment  of 

r as  given  by  (8),  the  variance  of  r can  be  found  to  be  given  by 


£[(l-a)2  + 2c2  + e2 


r2(l  -a)2 

K2 


(10) 


Comparing  (10)  with  (4.255),  we  note  that  for  2>a>0,  the 
convergence  of  the  mean  SNR  estimate  using  bias  uncorrected  esti- 
mates of  | T(F,tp) | ^ is  faster  than  if  these  biases  were  not 
corrected.  For  the  channels  considered  in  Section  4.1,  a.  is  a 
usually  positive  number  much  less  than  1 and,  therefore,  Eqs. 

(10)  and  (4.255)  are  essentially  equivalent. 

O 

In  summary,  the  use  of  unbiased  estimates  of  |T(F,t  )|  to 
estimate  the  mean  SNR  results  in  essentially  the  same  performance 
as  would  have  resulted  if  bias  corrected  estimates  of  |T(F,t  )|2 
were  used.  P 


A. 3 Estimation  of  the  RMS  Doppler  Spread 


The  rms  Doppler  spread  can  be  estimated  from  the  bias  uncor- 
rected estimates  of  jT(F,tp)|2  using  any  of  the  three  techniques 
analyzed  in  Section  4,3.2.  In  this  appendix  we  will  use  the 
differentiation  technique  to  determine  the  bias  in  estimating  B2 
due  to  using  bias  uncorrected  estimates  of  |T(F,t_)|2.  it  is 
felt  that  the  convergence  of  the  B2  estimator  will  not  be  appre- 
ciably affected  by  these  biases. 

2 

From  (4.169),  the  estimate  of  B for  the  differentiation 
technique  is  given  by 

/N  *T» 

2 

B - -o — h—  (H) 

it  (A tTT2 

where 

T1  “ K p^[l?(F>tp+A£)|2  ‘ |i(F>tp)|2]2 

T24^[lT(F,tp)|2]2  (12) 


Using  the  bias  uncorrected  estimates  of  the  magnitude 
squared  of  the  channel  transfer  function  as  given  in  (1),  then 
it  follows  that 


2 r 

E-jlJ  = 2|T(f,t)|2  (l-p0jAt)(l-c<)2  + 2<r2  + f2 


e{t2}  - 2|T(f,t)| 


2 

1 A.  & 

l-a  + — 
■» 


, 24.1S 

+ € +Tj 


4nJ  T|  c 4N 
+ -rVjr-  + 


VF>  P2(F) 

z 


(13) 


where  p is  defined  by  (4.156). 

The  mean  value  of  a can  be  determined  in  a manner  analogous 
to  that  used  in  Section  4.3.2.  Therefore,  for  K large  it  follows 
that 


A-4 


■.{?} 


0 2 , 2 
.y  0 2c  + €„ 

* a -a)  + -o | 

f x Atr 


1 4. a'  4.  2 €s  , 

1 - Qd  +*y  + e + “2“  + 


F + (f)/"2 


The  effect  of  not  correcting  the  bias  when  estimating 
|T(F,t D)i  is  t0  produce  a bias  in  estimating  even  when  the 
effects  of  c and  Co  are  negligible.  By  comparing  (14)  with 

p \ 

(4.193),  we  note  that  if  a and  — are  negigible  with  respect  to  1, 
then  bias  correction  may  not  improve  performance  significantly. 


A. 4 Estimation  of  RMS 


ipatt 


The  rms  multipath  spread  can  be  estimated  from  the  bias  un- 
corrected estimates  of  |T(F,tp)|2  using  any  one  of  the  three 
techniques  analyzed  in  Section  4.3.3.  For  the  differentiation 
techniques,  the  bias  in  estimating  L.2  is  the  dual  of  (14).  That 
is 


e{l2}  - 


2 2 

2 2 2< 

1/(1  -«r  + -5 — f 

V (&F ) 

l-*+4  +(2+T  + h(?f/2 


(15) 


Therefore,  if  a and  ^ are  negligible  with  respect  to  1,  then 
bias  correction  may  not  improve  performance  significantly. 


APPENDIX  B 


DETERMINATION  OF  BIAS  IN  ESTIMATING  THE 
BRANCH  ENVELOPE  CORRELATION  COEFFICIENT 


B.l  Introduction 

In  this  appendix  we  address  the  problem  of  determining  the 
bias  in  estimating  the  branch  envelope  correlation  coefficient. 
Four  cases  will  be  considered.  For  each  case  considered,  an  upper 
and  a lower  bound  on  the  expected  value  of  pjc<  will  be  found. 

The  expected  value  of  will  be  bounded  when  the  correlation 
coefficient  formula  (Figure  4.18)  or  the  simplified  formula 
[Eq.  (4.261)] is  used  with  or  without  correcting  for  the  biases 
in  es timating  |T(f,t)|^. 


B.2  Simplified  Approach  with  Bias  Correction 

From  (4.261)  of  Section  4.3,  the  branch  envelope  correlation 
coefficient  can  be  estimated  by 


T.  . 

* - , 

Pkj  ~ T.T.  “ 1 
i J 


(1) 


where 


lEji 


kj  Kp^L'k'^P 


K 


T,,(F,t„)|2(l+Ekp)  + 6kp  |TJ(F,tp)|2(l  + E,„)  + 6 


JP  JP 


\-kt  ivFVr<l+V  + #kP 

P 


1 A ? 

T.  - A y | T . (F , t ) | (1  + e . ) + 6. 

J K M 3 P JP  JP 


(2) 


The  bias  that  we  will  determine  is  the  asymptotic  bias;  that 
is , we  will  determine  the  bias  assuming  that  the  number  of  samples 
K is  sufficiently  large.  With  this  assumption,  the  mean  of  p.  . 
is  kJ 


B-l 


From  (2)  with  6 and  e uncorrelated, 


we  have 


E{\j}  ■ <1  + <Vj>'TJ(t*t»2  llTt(f.t)|2(l  + E{«kp«ip})+  KVjp} 

E{5kp8ip} 


E{W”ITj<f’t)|2  ITk(£>t)|2  1 + l'>kj(1+EKpfjP}) 


kp  ipJ 

K 

(4) 


For  large  K 


EjTjTj}  - ! Tj  (f  >t)i  |Tk(f,t)| 


(5) 


Since  we  can  bound  the  moments  of  the  errors  by 

0 s Ekp'JP} a <2 

0 s EkP5jp} £ fs  'TJc£*t>i2  ivf>t)|2  (6) 


Using  (4),  (5),  and  (6)  in  (3)  gives 


B . 3 Simplified  Approach  Without  Bias  Correction 

In  this  section  we  will  assume  that  the  biases  in  estimating 
|T(f,t)|2  have  not  been  corrected  (except  for  the  noise  bias;  see 
(2)  of  Appendix  A).  For  this  case,  (3)  is  still  valid  except 
that  the  required  expectations  are  given  by 


I 


E{TjTk}  “ 1Tj  C£,t)| 2 |Tk(f,t)|2 

E{\)}  - lTk(f.O|2  |Tj(f,t)l2[(l+pkJ)  l + E{tkpcjp}-2a+  2a(l  - a)J 

+ E{6kp6jp}  <8> 


Using  (6)  and  (8)  in  (3),  it  follows  that 


B.4  Formula  Approach  with  Bias  Correction 


In  this  section  we  will  use  the  correlation  coefficient 
formula  and  determine  the  bias  in  estimating  pk- . This  approach 
is  functionally  presented  in  Figure  4.18  of  Section  4.3.  The 


estimate  of  p^.  is  given  by 


Tki  - v- 


K-Tk)(vTj)r 


(10) 


where 


\4  eK<f>v,2(1 

p=lL  v 


Tk(F,tp)|2(l  + fkp)  +«kp' 


If  . ) + 6 . 
JP  JP 


(ID 


Using  bias  corrected  estimates  of  |T(F,t  )|Z  and  assuming  K 
is  sufficiently  large,  the  it  follows  that  (sle  Section  4.2): 
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__ 2 

K(Hk)  - E(Tk}  - i'ik(f,t)|2  (1.»-2c-2  + Cg) 


e{r,}  - e{Tj}  * “iTj  (f  ,t>|'2‘  (l  + 2c2  + eg) 


(12) 


From  (4),  (5),  (1.0),  and  (1.2)  we  have 


,Vj  + a+pkj)E{fkptjp}  + 


|Tk(f,t)h  | T . ( f , t*. ) ! 


2 2 
1 + 2e  + €g 


(13) 


Using  (6),  we  can  bound  E{$kj}  by 


Pu  -i 


1 + et  + 


rVrr  s E0Vj)  s 


2 (■: 


pki(i+f2)  + f2  + 4 

2 2 
1 + eg  + 2e 


(14) 


B , 5 Formu  la  Approach  thout:  bias  Corrcctioa 

When  the  formula  approach  given  by  (10)  is  applied  using 
bias  uncorrected  estimates  of  |T(ir,t)|^,  then  it  is  easy  to 
determine  that  for  K.  sufficiently  large 


'W  - 

e{r.}  - 


K)  “ lTk(£.o"l2  [d-“)2 

• Tj|  = |Tj  (f  ,t)|  2 [(1  - ot)2 


+ 


+ 


(15) 
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Using  (8)  and  (15),  the  expected  value  of  /5,  , is 


prjI1  + EKPejp}  - 2a)+  E{ekPfjp}  • 20,2 + 'k£ ^ 


Ekj} 


|T.(f,t)l2|Tk(f>t)l^ 


(1  - a)2  + 2c 2 + 


APPENDIX  C 


MINIMIZATION  OF  THE  RMS  ERROR  IN  ESTIMATING  |T(f,t)|2 
CONDITIONAL  UPON  THE  FADE  LEVEL 


In  this  appendix  we  address  the  problem  of  minimizing  the  rms 
fractional  error  in  estimating  the  magnitude  squared  of  the  chan- 
nel transfer  function  conditional  upon  a given  fade  level.  The 
utility  of  this  estimation  is  to  obtain  good  estimates  of  |T(f,t)|2 
at  a predetermined  fade  level.  From  (4.101)  of  Section  4.1.4  the 
variance  of  the  output  of  the  |T(f,t) |2  estimator  is  given  by 

i2  - f2|T00(F,t)|4  + f2|T(f,t)|22  (1) 


2 

where  |T(f,t) j has  not  been  normalized  to  1,  as  is  the  case  in 
(4.101). 

We  can  define  a fade  level  by 

,2 


A lToo^F,t^ 


F = ^ 

L | T(f , t)  | ^ 


In  Section  1 we  used  a fade  level  given  by  -10  log  F^< 


Using  (2)  in  (1),  it  follows  that 
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«2  - lT00(F.t)|4 


c + 2 

. FL 


(2) 


(3) 


The  rms  fractional  error  conditional  upon  a given  fade  level 
is  defined  as 


(4) 
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We  would  like  to  determine  the  filter  parameters  that  mini- 
mize (4).  Assuming  the  functional  form  of  the  filters  given  by 
(4.105)  and  with  it  follows  that 


./2  | 


+ 2L_. 
IT  WT  32 


+ -v/A, 


Comparing  (5)  with  (4.108)  of  Section  4.1.5,  we  note  that  the 
determination  of  W and  T such  that  (5)  is  minimized  is  identical  to 
the  minimization  of  (4.108).  In  fact,  it  should  be  noted  that  (5) 
and  (4,108)  are  identical  except  for  a scaling  of  B and  L. 
Therefore,  the  optimum  filter  parameters  can  be  found  from  (4.116) 
and  (4.117)  to  be  given  by 


rr  L2 


1 + A, 

(1  +A2^  TOl 


40 it  B L/ 
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/A1A2 
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(1+A1) 

1 
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L  1  2_ 

+ 1^A1A2  - !) 
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Comparing  (6)  and  (7)  with  (4.116)  and  (4.117)  we  note  that 
when  zq  is  minimized  for  a given  value  of  F^,  as  opposed  to  minimi* 

zing  t , W and  T change  by  a factor  Fj/3/2^6.  For  example, 
if  for^a  gi^/en  channll  it  is  desired  to  minimize  tr  for  FT  =0.001 


channel  it  is  desired  to  minimize  for  Fj.  =0.001 
(i.e,,  a 30  dB  fade),  then  W t and  T t are  approximately  one- 

tenth  of  the  values  that  minimized  e . 

P 
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Substituting  the  respective  optimum  filter  parameters  into 

(5)  and  (4.108),  then  eP  and  e can  be  related  by 

L p 


cC 


(8) 


In  summary,  optimum  filter  parameters  have  been  found  for 
estimating  |T(F,t)]2  at  a predetermined  fade  level.  Also,  (8) 
can  be  used  to  evaluate  the  rms  fractional  error  for  the  channels 
of  interest. 
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